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We show that there is exact dragging of the axis directions of local inerttal frames by a weighted 
average of the cosmological energy currents via gravitomagnetism for all linear perturbations of all 
Friedmann-Robertson- Walker (FRW) universes and of Einstein's static closed universe, and for all 
energy-momentum-stress tensors and in the presence of a cosmolgical constant. This includes FRW 
universes arbitrarily close to the Milne Universe and the de Sitter universe. Hence the postulate 
formulated by Ernst Mach about the physical cause for the time-evolution of inertial axes is shown 
to hold in General Relativity for linear perturbations of FRW universes. — The time-evolution of 
local inertial axes (relative to given local fiducial axes) is given experimentally by the precession 
angular velocity Hgyro of local gyroscopes, which in turn gives the operational definition of the 
gravitomagnetic field: Bg = — 2ngyro. The gravitomagnetic field is caused by energy currents via 
the momentum constraint, Einstein's G'\ equation, (—A 4- /i^)Ag = —IGttGnJs with Bg = curlAg. 
This equation is analogous to Ampere's law, but it holds for all time- dependent situations. A 
is the de Rham-Hodge Laplacian, and A = —curl curl for the vorticity sector in Riemannian 3- 
space. — In the solution for an open universe the 1/r^-force of Ampere is replaced by a Yukawa 
force Yp,{r) — (— d/dr)[(l/_R) exp(— /ir)], form-identical for FRW backgrounds with K = (—1,0). 
Here r is the measured geodesic distance from the gyroscope to the cosmological source, and 2nR is 
the measured circumference of the sphere centered at the gyroscope and going through the source 
point. The scale of the exponential cutoff is the T^-dot radius, where H is the Hubble rate, dot is the 
derivative with respect to cosmic time, and /.t^ = —4{dH/dt). Analogous results hold in closed FRW 
universes and in Einstein's closed static universe. — We list six fundamental tests for the principle 
formulated by Mach: all of them are explicitly fulfilled by our solutions. — We show that only 
energy currents in the toroidal vorticity sector with 1 = 1 can affect the precession of gyroscopes. 
We show that the harmonic decomposition of toroidal vorticity fields in terms of vector spherical 
harmonics A^ has radial functions which are form-identical for the 3-sphere, the hyperbolic 3-space, 
and Euclidean 3-space, and are form-identical with the spherical Bessel-, Neumann-, and Hankel 
functions. — The Appendix gives the de Rham-Hodge Laplacian on vorticity fields in Riemannian 3- 
spaces by equations connecting the calculus of differential forms with the curl notation. We also give 
the derivation of the Weitzenbock formula for the difference between the de Rham-Hodge Laplacian 
A and the "rough" Laplacian on vector fields. 

PACS numbers: 04.20.-q, 04.25.-g, 98.80.Jk 



I. SUMMARY AND CONCLUSIOS 

A. Conclusions 

The time-evolution of local inertial axes, i.e. the local 
non-rotating frame, is experimentally determined by the 
spin axes of gyroscopes, e.g. by inertial guidance systems. 
This has been first demonstrated by Leon Foucault in 
1852. 

It is an observational fact that the spin axes of gyro- 
scopes far away from Earth do not precess relative to 
quasars. This observational fact has been named "Mach 
zero" by Bondi and Samuel [l|. — Near Earth there are 
two corrections predicted by General Relativity: (1) The 
de Sitter precession, which is due to the gyroscope's mo- 
tion in the curved Schwarzschild metric, and which has 
been measured. (2) The extremely small frame-dragging 
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effect by the rotating Earth, the Lense-Thirring effect, 
which is predicted to be 43 milli-arc-sec per year, which 
hopefully will be extracted from the data taken by Grav- 
ity Probe B 2], and whose measurement by the LAGEOS 
satellites has been reported in Q . — The dragging effect 
by the rotating Earth (and also the perihelion shift of 
Mercury) is measured relative to distant quasars. There- 
fore a measurement by Gravity Probe B and the LA- 
GEOS satellites (and a measurement of the perihelion 
shift) is a test of two things combined, on the one hand a 
test of Einstein's General Relativity, on the other hand 
a test of the observational fact "Mach zero" . 

But the fundamental question remains: What physi- 
cal cause explains the observed time-evolution of gyro- 
scope axes, the observational fact "Mach zero"? In J. A. 
Wheeler's words: Who gives the marching orders to the 
spin axes of gyroscopes, i.e. to inertial axes? Can theory 
predict resp. explain the extremely precise observational 
tests of "Mach zero"? 

An answer to this fundamental question was formu- 
lated by Ernst Mach 0, in his postulate that inertial 
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axes exactly follow an average of the motion of masses 
in the Universe. Mach's postulate demands exact frame 
dragging, not merely a little bit of frame dragging as in 
the Lense-Thirring effect. — Mach did not have a mech- 
anism for dragging. The new force came with General 
Relativity: Gravitomagnetism, exemplified in the Lense- 
Thirring effect. — Mach also did not know, what average 
of the cosmological mass currents should be taken. 

Mach's starting postulate was: No absolute space, in 
modern parlance no absolute element in the input. Only 
relative motion of physical objects is significant. — From 
this requirement follows Mach's postulate of frame in- 
variance of observable predictions under time-dependent 
rotations of the reference frame. In Mach's words 0: 

"Obviously it does not matter, whether we 
think of the Earth rotating around its axis, or 
if we imagine a static Earth with the celestial 
bodies rotating around it." 

From these requirements follows Mach's postulate j^: 
Centrifugal forces and the rotation of the plane of Fou- 
cault's pendulum are totally determined by the relative 
motion of masses in the universe with respect to Earth- 
fixed axes. In modern parlance: The time-evolution of 
local inertial axes is totally determined by the energy cur- 
rents in the universe including the effect of gravitational 
waves. — The demands "no absolute element," "frame 
invariance under rotational motion," and "totally deter- 
mined" necessitate the postulate of exact dragging of in- 
ertial axes by energy currents in the universe. With only 
partial dragging it would be impossible to fulfill any of 
the three postulates " no absolute element in the input" , 
"frame invariance of the solutions", and "time evolution 
of inertial axes totally determined by cosmological energy 
currents" . 

For a purely local analysis, if one has only laboratory 
experiments available (like Newton's bucket experiment), 
the local nonrotating frame is an absolute element. — 
But, in contrast, for a global analysis, if one considers 
observations of the entire universe, the local nonrotating 
frame is not an absolute element in the input, because 
(according to Mach) the nonrotating frame is totally de- 
termined by mass motions (energy currents) in the uni- 
verse, it is nonrotating relative to an average of the energy 
currents in the universe. 

The "relativists" (Huyghens, Berkeley, Leibniz, Mach, 
and others) wanted a theory without absolute space 
(without an absolute element in the input), a theory, 
where only relative motion is significant. After Ein- 
stein's theory of General Relativity it was soon recog- 
nized that this theory, when applied to the solar system, 
is not a "Theory of Relativity" in the sense that the solu- 
tions violate the postulates "no absolute element" , "solu- 
tions frame invariant" , "totally determined" , and "exact 
dragging" , see Sec. II HI This is why Einstein and oth- 
ers started to focus on cosmological models. The un- 
perturbed Friedmann-Robertson- Walker models satisfy 
Mach's postulates. In this paper we derive the new result 



that all FRW models with arbitrary linear perturbations 
also satisfy Mach's postulates. This means that Cosmo- 
logical General Relativity (at the level of linear perturba- 
tions of FRW universes) is a true "Theory of Relativity". 

We have shown in [y, 0] that the axis directions of 
local inertial frames everywhere in the universe exactly 
follow a weighted average of the energy currents in the 
universe, i.e. there is exact frame- dragging, for linear 
perturbations of a Friedmann-Robertson- Walker (FRW) 
background with K = (i.e. spatially flat). The weight 
function has an exponential cutoff at the H-dot radius, 
where H is the Hubble rate, and dot is the derivative 
with respect to measured cosmic time. Hence we have 
demonstrated the validity of the hypothesis formulated 
by Mach about the physical cause for the time-evolution 
of inertial-frame axes for linear perturbations of a FRW 
background with K ^ 0. The proof follows from Ein- 
stein's Gqj equation, the momentum constraint. The dy- 
namical mechanism is cosmological gravitomagnetism. 

In the present paper we show that results analogous 
to those derived for X = in [3] are valid for linear 
perturbations of FRW backgrounds with K = ±1, i.e. 
backgrounds which are spatially spherical (S^) resp. hy- 
perbolic (H^). 

Many alternative and inequivalent postulates have 
been proposed under the heading "Mach's principle" by 
many authors after E. Mach. An incomplete list of 10 in- 
equivalent versions of Mach's principle has been given by 
Bondi and Samuel [Ij] (based in part on fsj). While the 
versions "Mach 9" and "Mach 10" listed by 1] are closely 
related to our "Fundamental Tests 3 and 4" in Sec. IV HI 
their versions "Mach 1-8" definitely do not correspond 
to anything in the writings of Mach. Furthermore each 
one of their versions "Mach 1-8" is either not true in cos- 
mological General Relativity, as stated already in or 
irrelevant, as we explain in Sec. II Jl 

A more detailed introduction to the conceptual issues 
of the postulate formulated by Mach and discussions of 
publications @-[l^ (niostly cosmological as postulated 
by Mach) are given in Sees. I and II of the companion 
paper 0. Additional references are given in 8, 12, ij, 
IitI . ITsj . Comparisons of our methods and results with 
the recent papers by Bicak, Lynden-Bell, and Katz [l^, 
[III are at the end of Sects. |TT] and IVTll — We use the 
conventions of Misner, Thorne, and Wheeler [20j . 

B. Vorticity perturbations 

We now summarize the main results of the present pa- 
per for the case when only vorticity perturbations (= 
3-vector perturbations, by definition divergenceless) are 
present. 

Two important results from cosmological perturbation 
theory [21] for K = (0, ±1) are needed to understand the 
following summary. In the vorticity sector: 

1. The slicing of space-time in slices Ei of fixed time, 
i.e. 3-spaces, is unique. The lapse function (elapsed 
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measured time between slices) and goo are unper- 
turbed. 

2. The intrinsic geometry of each slice Et, i.e. of 3- 
space, remains unperturbed. 

Details and proofs are given in Sec. Ill of 7]. Note that 
Sees. Ill - V of are written for FRW backgrounds with 
K={0,±1). 

For the 3-geometry, which remains unperturbed {E^, 
resp. S^, resp. H^), and for our fixed-time problem (the 
momentum constraint), the uniquely appropriate coor- 
dinate choice (gauge choice) is comoving spherical coor- 
dinates, = (x,0,</)), where (ax) is the measured ra- 
dial distance with a{t) = scale factor. With this gauge 
choice the coefficients of the 3-metric gij are also un- 
perturbed. The spatial coordinate system is geometrically 
rigid apart from the uniform Hubble expansion. Also goo 
is unperturbed, and we choose the time-coordinate t to be 
the measured cosmic time. Hence, with our gauge choice 
the only quantity referring to vorticity perturbations is 
the shift 3- vector /?' (resp Pi = goi), 

ds^ = ~dt^ + a^[dx^ + Rlomid9^ +s^'^^&d<j)^)] 

+2(3^ dx'dt, (1) 

where i?com stands for i?comoving, and i?com(x) = 
(x, sinx, sinhx) for K = (0, +1,— 1). The index of the 
shift 3-vector is lowered and raised by the 3-metric, e.g. 

Our choice of unperturbed comoving spherical coordi- 
nates for {E^, S'^, H^) is singled out by the fact that our 
radial coordinate lines {9 and 4> fixed) are geodesies in 
3-space, an important property in discussions of Mach's 
principle. In all other gauges, where perturbed 3-metrics 
g^v are used, the radial coordinate lines wind up more 
and more into "spirals" (as time goes on) relative to 
geodesies on S^. Hence all other gauges (particularly 
time-orthogonal) give a very awkward way to coordina- 
tize unperturbed {E^,S^,H^) spaces. 

Up to now the coordinatization (gauge) is not yet com- 
pletely fixed, there still exists the freedom of residual 
gauge transformations by time-dependent, spatially rigid 
rotations of the coordinates for 3-space. Rigid rotations 
are crucial in Mach's statement on frame- invariance [![: 
See the quote in Sec. lA, "Obviously it does not matter 

We shall show in Sec. II HI that our solution of the equa- 
tions of cosmological gravitomagnetism is explicitly form- 
invariant under time-dependent, spatially rigid rotations 
of the reference frame, i.e. frame-independent (valid in 
all reference frames connected by time-dependent rigid 
rotations). In contrast, the solutions of the equations 
of General Relativity for the solar system are not form- 
invariant under spatially rigid, time-dependent rotations. 

Measuring the angular velocity of a given celestial body 
requires that one first fixes the residual gauge freedom, 
i.e. one must first state, relative to what one wants to 
measure the angular velocity. This can be done by mak- 
ing a definite choice of spatial axis directions along one 



world line. In the spirit of Mach one could fix the axis di- 
rections either along the world-line of the Earth to Earth- 
fixed axes, or (in an open, asymptotically unperturbed 
universe) to asymptotic quasars. 



C. Gravitomagnetic and gravitoelectric fields 

The general operational definitions (valid beyond per- 
turbation theory) of the gravitomagnetic field Bg and 
the gravitoelectric field Eg are given via measurements 
by fiducial observers (FIDOs) with nonintersecting world 
lines and with their local ortho-normal bases, LONBs 
[2^ . These FIDOs measure the first time-derivatives of 
the LONB components of the 3-momentum, , for free- 
falling quasistatic test particles, resp. of the spin-vector. 
Si, for comoving gyroscopes, 

free-falling quasistatic test particle: — (pj) = mEf, (2) 
gyro comoving with FIDO: ^(^.) = -\{Bg x S)^, (3) 

I 2 ^ 

where Q,^ are the LONB components of the precession 
rate of the spin axes of a gyroscope. Hats over indices de- 
note components in a local ortho- normal basis (LONB), 
and t is the local time measured by the FIDO. Note that 
LONB components, e.g. pj, are directly measurable (for 
a given field of LONBs), while coordinate-basis compo- 
nents, e.g. p^, are not measurable without prior knowl- 
edge of the metric g^,^. 

The operational definitions of Ba and Ea are identical 
with the operational definitions of the ordinary electric 
and magnetic fields via measurements on charged point 
test-particles and on charged spinning test particles, ex- 
cept that the charge q is replaced by the mass m for 
quasistatic test particles. 

The (Bg, i5g)-fields depend on the choice of FIDOs. 
For free-falling FIDOs Eg = 0, and for nonrotating FI- 
DOs (relative to axes of local gyroscopes) Bg = 0, which 
are the twin conditions in the equivalence principle. — 
The {Eg,Bg) fields are directly related to connection 
coeflacients: (^^5)0 = -Ef, (w-)6 = ^l^'f ^^ere 
Bvj = ^IjkBy ^Ij = ^Ijk^y and = (wj^)^. More de- 
tails about the general operational definitions of {Eg, Bg) 
are given in Sec. IV of [7|]. 

Our specific choice for a field of FIDOs: Our FIDOs 
are at fixed values of our coordinates x*, and we fix the 
spatial axis directions of our FIDOs in the directions of 
our coordinate basis vectors {d/dx, d/d9, d/dcf)). — We 
recaU from Sec. II Bl that our choice of unperturbed co- 
moving spherical coordinates for {E^,S^,H^) is singled 
out by the fact that the radial coordinate lines {9 and 
fixed) are geodesies in 3-space. Additional details about 
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our specific choice of FIDOs are given in ^ at the end 
of Sec. IV. 

The gravitomagnetic vector potential Ag in the vortic- 
ity sector is uniquely determined by Bg = curing, be- 
cause div Ag = 0. — With Pi = goi and Af defined 
independently above, it follows that the shift vector [3 is 
equal to the gravitomagnetic vector potential Ag, see also 
Sec. V of — Our gravitomagnetic vector potential 
Ag is directly proportional to Bardeen's gauge- invariant 
amplitude ^' as shown in Eq. (|13p . 

In Sec. mil we shall consider gravitomagnetism on a 
Minkowski background and show that our FIDOs (and 
hence our Eg and Bg fields) are singled out by the prop- 
erty that all 13 nonzero components of the Riemann ten- 
sor are uniquely determined by spatial and time deriva- 
tives of our Bg field alone. Our Bg field gives an efficient 
and physically transparent way to represent all informa- 
tion in the Riemann tensor. — In the opposite direction: 
We shall show that our Bg field is determined by the Ricci 
components Rj^Q up to a homogeneous Bg field, which is 
equivalent to a time-dependent rigid rotation of the co- 
ordinates and FIDO axes according to Eq. 

D. The momentum constraint 

The crucial equation for Mach's principle is the mo- 
mentum constraint, Einstein's equation, for cosmo- 
logical vorticity perturbations. The momentum con- 
straint for perturbations of all three FRW background 
geometries, K — (±1,0), is given by one form-identical 
equation without curvature terms: 

(-A + /?)Ag = -IGttGn Je, (5) 

where ifJ-/2)^ = -(dH/dt) = (Hubble-dot radius)-2 
Derivation of Eq. ^ in Sec. [Til Our momentum con- 
straint has the form of the inhomogeneous Helmholtz 
equation with the sign of the /i^-term opposite to the 
usual sign. The sign in Eq. (O causes an exponential 
cutoff, which becomes important beyound the Hubble- 
dot radius. — The source with Jf = for all types of 
matter is the directly measurable energy-current density, 
which is equal to the measurable momentum density (for 
c = 1), hence the name "momentum constraint". "Di- 
rectly measurable" means that J~ is measurable without 
prior knowledge of goi; which is the output of solving 
the momentum constraint. — The energy current is 
analogous to the charge current Jq in Ampere's law. For 
perfect fluids and linear perturbations — {p -\- p)v. — 
All symbols in Eq. ([5]) refer to the physical scale, they 
are not rescaled to a comoving length scale. Therefore 
the scale factor a{t) does not appear in Eq. ([5]). 

It is remarkable that Eq. ([5]) holds for time- dependent 
gravitomagnetodynamics. All the same, it is an elliptic 
partial differential equation: No partial time-derivatives 
of the perturbation field Ag. Time derivatives only appear 



in the given background, H = a ^da/dt, H = dH/dt. — 
Eq. ([5]) has the same form as Ampere's law for stationary 
magnetism except for the p.^Ag term. — The analogue of 
the Maxwell term {dtE) in the Ampere-Maxwell equation 
is absent for the time- dependent context of gravitomag- 
netism. No inconsistency arises when taking the diver- 
gence of Eq. ([5]), since div Ag = and div = in the 
vorticity sector. The analogue of a Maxwell term can- 
not be present in gravitomagnetism, this would produce 
gravitational vector waves, which is impossible. Addi- 
tional explanations about the momentum constraint are 
given after Eq. (36) in Sec. VI of [7]. 

Other perturbed Einstein equations and the perturbed 
evolution equations for matter will not be needed in the 
context of this paper. The full set of field equations for 
linear cosmological gravitomagnetism is given in Sees. V 
and VI of 7]. 



E. The de Rham-Hodge Laplacian on vector fields 
in Riemannian 3-spaces 

The absence of a 3-curvature term in Eq. ((S]) is con- 
nected to our use of the de Rham-Hodge Laplace oper- 
ator A, which is called "the" Laplace operator in the 
literature on differential forms in Riemannian geometry 
[l^-jl^. Except when acting on scalar fields, the de 
Rham-Hodge Laplace operator A must be distinguished 
from = S^^V^V^, named the "rough Laplacian" in 
[H, HEl- — Unfortunately has been used in all publi- 
cations on cosmological vector perturbations (as far as we 
know), e.g. in [l!], [23], and moreover has been called 
Laplacian and/or denoted by A in many publications on 
vector perturbations, e.g. in [2^-[3l|. Using has the 
unfortunate consequence that curvature terms appear, 
where none are present when the de Rham-Hodge Lapla- 
cian A is used. The most important example is electro- 
magnetism in curved space-time, where the equivalence 
principle forbids curvature terms, if the de Rham-Hodge 
Laplacian resp d' Alembertian is used [111 . 

In Riemannian 3-space, the de Rham-Hodge Laplacian 
acting on vorticity fields, i.e. fields with div a = *d-kd = 
6a = 0, is given by 

(Aa)^ — —(curl curl a)p ~ —(-kdi^da))^, (6) 

where we have given both the notation of elementary 
vector calculus and the notation of differential forms with 
d = exterior derivative, {—d) = d* = adjoint operator to 
the exterior derivative, ★ = Hodge dual, tilde ^ denoting 
p-forms, and (a)^ = (a)^. Note that (curla)^ = (-kdd)^. 
Details are given in the Appendix „ 

The de Rham-Hodge Laplacian is singled out among 
the elliptic differential operators by the following prop- 
erties: (1) If in some region all types of sources are zero, 
div V — and curl v = 0, then also A v = 0. (2) 
A commutes with d, 6 and Therefore applying curl, 
div, and grad to an eigenfield of the Laplacian produces 
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again an eigenfield of the Laplacian with the same eigen- 
value. — The "rough Laplacian" does not have prop- 
erties (1) and (2), if the Ricci tensor is not zero. — The 
de Rham-Hodge Laplacian is also easier to compute than 
the "rough Laplacian", because no Christoffel symbols 
resp. connection coefficients are needed to compute it. 

The difference between the de Rham-Hodge Laplacian 
and the "rough Laplacian" is given by the Weitzenbock 
formula, i.e. by curvature terms. The Weitzenbock for- 
mula for vector fields, is derived in Sec. [21 of the Ap- 
pendix. For a FRW universe with K — ±1 and with the 
curvature radius Cc = 1, the Weitzenbock formula gives 
AA = (V2 - 2K)A. 

The de Rham-Hodge Laplacian applies to totally an- 
tisymmetric tensor fields. For totally symmetric ten- 
sor fields, e.g. for cosmological tensor perturbations, 
the appropriate Laplacian is the Lichnerowicz Laplacian 
[33I [3^ . For vector fields the Laplacians of Lichnerowicz 
and of de Rham-Hodge coincide. 



F. No absolute element in the input: 
The exponential cutoff at the //-dot radius 

The crucial difference between our momentum con- 
straint, Eq. ([5|), and the corresponding equation of 
Bardeen, Eq. (4.12), and the equation advocated by 
Bicak et al [H, [l^ is our fi^Ag term, which causes the 
Yukawa cutoff in Eq. ^ below, and which is absent in 
the solution of the equation advocated by by Bicak et al. 
Our n^Ag term is due to the difference of the sources. 

Our guiding principle and starting point: "No Ab- 
solute Element in the Input". The source, the matter 
input, must be measurable without prior knowledge of 
goi = aI^"* . The latter is the output of solving the momen- 
tum constraint. From Eq. Q we see that Bg = curlAg 
directly gives the local nonrotating frame all over the 
universe, which should not be in the input, because it 
would be an absolute element in the input. — We write 
the momentum constraint in terms of the LONE compo- 
nents T*?, which are measurable without prior knowledge 

of Ag, i.e. without an absolute element in the input. 

For Mach's principle, the relevant factor in the source 
is the angular velocity of cosmic matter, rimatter, mea- 
sured from our position and relative to geodesies (on T,t) 
which start in the directions of our chosen fiducial axes. 
In the local ortho-normal bases, the transverse part of 
the energy current ~ ji"'' is directly proportional 
to £7™''**°'' (relative to our fiducial axes), which can be 

measured without prior knowledge of goi = A] . The 
measured transverse matter velocity is directly propor- 
tional to Bardeen's gauge-invariant velocity amplitude 
Vs, which he does not use in his momentum constraint. 

In contrast, Bardeen's Eq. (4.12) (and the equation ad- 
vocated by Bicak et al) is written in the coordinate basis 
with the source T^, and this is proportional to Bardeen's 



other gauge- invariant velocity amplitude Vc- The curl of 
the velocity field with the amplitude Vc is proportional to 
the angular velocity of matter relative to the axes of gy- 
roscopes at the positions of the cosmic sources. But this 
has the severe drawback that Vc is not measurable without 
prior knowledge of the time-evolution of gyroscope axes 
all over the cosmos, which should be the output of solving 
the momentum constraint: r^gyro ~ ^5^g = ~'^^^^^^g 
with A^f^ — goi. Bardeen's version of the momentum 
constraint (and the version advocated by Bicak et al) 
needs an absolute element in the input: the local non- 
rotating frames all over the cosmos. This violates the 
crucial starting demand: No absolute clement in the in- 
put. The requirement "the theory should not contain an 
absolute element" , discussed in [33 , is closely connected 
to Mach's starting point: No absolute space. 

Einstein commented about the energy-momentum ten- 
sor T^i, in a coordinate basis: "If you have a tensor T^^ 
and not a metric ... the statement that matter by itself 
determines the metric is meaningless" [s^. This objec- 
tion of Einstein directly applies to the approach advo- 
cated by Bicak et al nq . . On the other hand we have 
shown in Sec. IX of Q that this objection of Einstein 
does not apply to the LONE components of the energy- 
momentum tensor. 



G. The solution of the momentum constraint 

The solution of the momentum constraint Eq. ([5]) for 
cosmological gravitomagnetism on a background of open 
FRW universes is given by identical expressions for K = 
and K = —1, 

Bg{P) - -4Gn j d[Y0\Q)[npQ X UQW^.{rpQ), (7) 

Yuir) = — [ — exp(— /ir)] = Yukawa force, (8) 
dr R 

where rpq is the measured geodesic (radial) distance 
from the field point P (gyroscope) to the source point Q, 
and [211 Rpq) is the measured circumference of the great 
circle through Q with the center P. Hence r = ax, where 
a is the scale factor, and x is the radial coordinate (co- 
moving geodesic distance from the origin). — The struc- 
ture of Eq. ([7]) is analogous to Ampere 's force for ordinary 
stationary magnetism, however the of Ampere's 

force is replaced by the Yukawa force of Eq. ([8]). — In 
the curved space the geodesic distance rpq appears 
in the exponential cutoff and in the derivative of Eq. ([5]), 
while RpQ appears in the denominator. 

Although the structure of Eq. ([7]) is analogous to the 
solution of Ampere's law for ordinary stationary mag- 
netism, there is a fundamental difference to Ampere's 
law, since Eq. ([7]) is valid for time- dependent situations, 
i.e. gravitomagneto-dynamics. 

The vector structure of Eq. ([7]) needs explanation, since 
there is no global parallelism on a hyperbolic 3-space . 
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The vectors Je{Q) and npq are in the tangent space at 
the source point Q, and npQ is the unit vector pointing 
along the geodesic from P to Q. We parallel transport 
[npQ X Jj(Q)] along the geodesic from Q to P, where the 
contributions from all sources are added. — We use the 
arrow-notation of Thorne et al [l^l for vectors in tangent 
spaces to Riemannian 3-spaces from a (3 + l)-split. This 
emphasizes the structure analogous to Ampere's law. 



H. Frame invariance, no boundary conditions 

There is a second fundamental difference between the 
solution Eq. ([7]) for cosmological gravitomagnetism and 
the corresponding solutions in other theories, Ampere's 
magnetism, electromagnetism in Minkowski space, and 
general relativity for the solar system. The solutions 
in all these other theories contain an absolute ele- 
ment: Ampere's law and the laws of electromagnetism 
in Minkowski space do not hold in frames which are ro- 
tating relative to inertial frames, unless one introduces 
fictitious forces. 

Einstein's equations by themselves are form- invariant. 
But for general relativity of the solar system, when work- 
ing in rotating frames relative to the asymptotic nonro- 
tating frame (the absolute element), one must impose 
boundary conditions on the solutions in the asymptotic 
Minkowski space. The boundary conditions encode the 
absolute element, they explicitely encode fictitious forces 
(which do not have sources within the solar system), 
as explained in Sec. XI of 01 • Hence the solution of 
Einstein's equations for the solar system is not form- 
invariant when going to a frame which is rotating relative 
to asymptotic inertial frames. 

In contrast, the solution for cosmological General Rel- 
ativity, Eq. ([7]), remains valid as it stands, i.e. the solu- 
tion is form-invariant, if one goes to a frame which is in 
globally rigid rotation relative to the previous reference 
frame: The form of the solution is frame-independent, 
because both sides of Eq. ([7]) change by the same term. 
Therefore the solution Eq. ([7]) contains no absolute ele- 
ment. — No boundary conditions at spatial infinity are 
needed for regular solutions in Eq. ([7]). 



From Eqs. (|4l [7]) we obtain 

/>oo 

f^gyro = < flrr^er >^ dr (i'ZZAr) W{r),{Q) 

Jo 

Wir) = ^l6TrGn{p + p)R''Y^{r) (10) 

for perturbations of open FRW universes. Y^{r) is the 
Yukawa force given in Eq. ([8]). 

The crucial equation for establishing exact (as opposed 
to partial) dragging of local inertial axis directions is the 
condition that the weight function W{r) must be normal- 
ized to unity, 

/>oo 

/ drW{r) = l, (11) 

as it must be for a proper averaging weight function in 
any problem. Our weight function with its Yukawa force, 
Eq. Uni), fulfills Eq. (HH) as shown in Sec. IVCl 

Eqs. ([HI fTT|l state our most important result: figym is 
equal to the weighted average (with proper normalization 
to unity) of rimattcr('') for all types of matter-energy and 
for all types of energy-current distributions. The time- 
evolution of inertial axes exactly follows the weighted av- 
erage of cosmic matter motion. Hence the evolution of in- 
ertial axes is fully determined (not merely influenced) by 
the cosmic energy currents. There is exact dragging (not 
merely partial dragging) of inertial axes by the weighted 
average of cosmic energy currents, as postulated by Ernst 
Mach Mach had asked: 

"What share has every mass in the determi- 
nation of direction ... in the law of inertia? 
No definite answer can be given by our expe- 
riences." 

Exact dragging of inertial axes by arbitrary cosmologi- 
cal energy currents in linear perturbation theory has been 
demonstrated for the first time in our papers 0, 0| for 
K = and in the present paper for X = ±1. 

J. Alternative, inequivalent postulates proposed by 
authors after Mach 



I. Exact dragging of inertial axes 

The solution Eq. ([7]) can be rewritten to show exact 
dragging of inertial axes explicitly. For reasons of sym- 
metry under rotations and space reflection relative to a 
gyroscope at the origin, a general velocity field of mat- 
ter for rpQ fixed can contribute to the gyroscope's pre- 
cession only through its term with (^ = 1, P = -I-) in 
the spherical harmonic decomposition. This term is a 
toroidal vorticity field, and it is equivalent to a rigid ro- 
tation with the angular velocity ri^^j^j,j.(r), see Sec. lV Al 



Our presentation of Mach's principle closely agrees 
with the one of Misner, Thorne, and Wheeler [s^l and 
Weinberg [3l|. However, many alternative and inequiv- 
alent postulates have been proposed under the heading 
"Mach's principle" by many authors after E. Mach. An 
incomplete list of 10 inequivalent versions of Mach's prin- 
g3le has been given by Bondi and Samuel [l|, see also 

In the list of Bondi and Samuel the version "Mach 9" 
( "The theory contains no absolute elements" ) is the pos- 
tulate of the so-called "relativists": Huyghens, Bishop 
Berkeley, Leibniz, Mach, and others. It was the starting 
point for Mach. See our "Test 3" in Sec. IVHl — The 
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version "Mach 10" ("Overall rigid rotations of a system 
are unobservable" ) is a necessary condition, formulated 
by Mach in his earliest writings Q quoted in Sec. II Al 
See our "Test 4" in Sec. IV HI Globally rigid rotations 
are discussed in Sec. IV Gl 

In contrast the versions "Mach 1, 2, 6, 7" are not true 
in cosmological General Relativity, as already stated in 
[l[. — The version "Mach 3" ("local inertial frames in- 
fluenced ...") is irrelevant, because with only partial 
dragging one violates Mach's postulates "no absolute el- 
ement" in the input and "frame-invariance" in the solu- 
tions (see the quote of Mach in Sec. LA). — "Mach 4" 
("The Universe is spatially closed") is irrelevant, because 
we show that exact dragging also holds in open uni- 
verses. — "Mach 5" ("The total angular momentum of 
the universe is zero" ) is not true for that angular momen- 
tum which can be defined and is observable without an 
absolute element in the input (i.e. without prior knowl- 
edge of the local non-rotating frames), see Sec. lIFI and the 
end of Secini — "Mach 8" ("p/pcrit is of order unity") is 
irrelevant, because we show in Sees. IV Dl IV Fl that exact 
dragging also holds for p/pcrit ^ 1 and p/pcrit *C 1- 

With respect to this list of 10 versions of Mach's prin- 
ciple we conclude: The criticisms "There are many for- 
mulations of Mach's principle" and "It is not clear, what 
Mach meant" do not apply. 

K. Outline 

In Sec. iniwe derive the momentum constraint for (1) 
that momentum density which can be measured as an 
input without prior knowledge of the solution of the mo- 
mentum constraint (the output), and for (2) the Lapla- 
cian A of de Rham-Hodge as opposed to the "rough 
Laplacian" V^. 

In Sec. mil we show that our Bg field on a Minkowski 
background uniquely determines all nonzero components 
of the Riemann tensor. 

To derive the Green function for the momentum con- 
straint for the curved 3-spaces and H^, we need the so- 
lution in source-free regions, i.e. we need vorticity eigen- 
flelds of the Laplacian. 

In a first step, in Sec. lIV A] we derive scalar eigenfunc- 
tions of the Laplacian on and H^. Our new result: We 
give the radial functions of the scalar harmonics in a sim- 
ple, identical form for all three geometries {E^ , , H^). 

In Sees. IIV Bl - IIV D] we give the expansion of vorticity 
fields (toroidal and poloidal) in terms of vector spheri- 
cal harmonics and Yirn^x- We show that the radial 
eigenf unctions of the Laplacian for toroidal vector fields 
are the spherical Bessel, Neumann, and Hankel functions 
generalized to and by the simple replacements 
r (x resp R) given in Sec. II Gl 

In Sec. |V] we show that only toroidal Ag-fields with 
^ = 1 are relevant for the precession of gyroscopes at the 
origin and that the precession of a gyroscope cannot be 
caused by scalar or tensor perturbations. We then derive 



the radial Green functions for the operator (A — /z^) act- 
ing on toroidal Ag-fields with i = I on {E'^ , ,H^). The 
final result is the gravitomagnetic Bg-field at the position 
of the gyroscope, i.e. the precession of the gyroscope, 
expressed by the energy currents in the universe, and 
the equation which expresses exact dragging of inertial 
axes. We also show that Mach's principle remains valid 
for perturbations of a universe with arbitrarily small to- 
tal energy inside the Hubble radius, i.e. arbitrarily close 
to the Milne universe, and for a universe arbitrarily close 
to a de Sitter universe. — We derive the gravitomag- 
netic fields for the zero modes of the operator (A -I- 4/C), 
i.e. for time-dependent globally rigid rotations of coordi- 
nates and FIDO axes in a FRW universe. — In Sec. IV HI 
we give Six Fundamental Tests for Mach 's Principle, all 
intimately related to Mach's original formulation and all 
explicitely fulfilled by our solutions. 

The Appendix n gives the tools for the de Rham - Hodge 
Laplacian on vector fields in Riemannian 3-spaces, specif- 
ically the explicit equations connecting the calculus of 
differential forms with the elementary notation with curl 
etc. These tools are essential for cosmological gravito- 
magnetism. 

II. THE MOMENTUM CONSTRAINT 

The momentum constraint written in local orthonor- 
mal tetrad components is given in Eq. ([5]) . It can be de- 
rived via two different general methods: (i) One can use 
the local orthonormal tetrad (LONB) method of Cartan 
from beginning to end. This method has been presented 
in Sees. V and VI of [Tj], where the connection coeffi- 
cients for vorticity perturbations are given for FRW back- 
grounds with K = (0, ±1), and Ggj is derived for K = 0. 
We have used this method to derive Eq. (O also for 
K — ±1, but we shall not present details again, (ii) Al- 
ternatively one can, in a first step, take the detour via the 
standard coordinate-component method with Christoffel 
symbols to obtain the momentum constraint in a coor- 
dinate basis with as the source. R^^ has been de- 
rived in _Bardeen [2l|, Eq. (A2b), and in Kodama and 
Sazaki |2^, Eq. (D.14b). In a second step one applies the 
basis transformation from the coordinate basis to the lo- 
cal orthonormal basis. — A modification of method (ii) is 
quicker: The basis transformation at the end of method 
(ii) can be replaced by taking Bardeen's momentum con- 
straint and moving the difference (T" — T"), which is 
meaningful at the level of linear perturbation theory, to 
the left-hand side of the momentum constraint. 

The momentum constraint for i?° of Bardeen and of 
Kodama and Sazaki reads in our notation 

(V2 + 2K/al)Al = 167rGN(po + PoK, (12) 

where Gc denotes the curvature scale, and refers to the 
physical scale. The 3- vector index i in and w* is low- 
ered and raised by the 3- metric. — The zero modes of the 
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operator {A+AK / a^) correspond to time-dependent rigid 
rotations of the coordinates as demonstrated in Sec. IV Gl 
For non-zero modes of (A + AK / ac) , our gravitomagnetic 
vector potential Ag, which is equal to our shift vector field 
is directly proportional to Bardeen's gauge- invariant 
amplitude ^ (a combination of shift amplitude and shear 
amplitude) , 

Alix, t) = P\x, t) = -vl/(t)Q(i)Xa;), (13) 

for a given wave number and polarization. For zero 
modes there is no shear, and a nonvanishing gauge- 
invariant amplitude does not exist. — Similarly for 
Bardeen's gauge-invariant velocity amplitude Vc 

vl{x,t)^v,{t)Q^'^\x). (14) 

To convert Bardeen's momentum constraint Eq. IT^ 
to our momentum constraint, the first step is replacing 
the "rough Laplacian" — V'Vi by "the" Laplacian 
A (de Rham-Hodge Laplacian) via the Weitzenbock for- 
mula. With the Ricci tensor of the 3-space for FRW the 
Weitzenbock formula gives Eq. (jA.19[) . 

{A-V^)A^-{2K/al)A. (15) 

Hence Bardeen's momentum constraint written in terms 
of the de Rham-Hodge Laplacian is 

(A -I- 4K/al)Al = WnG^ipo + PoH- (16) 

The second step is replacing the gauge-invariant mat- 
ter velocity amplitude ?;c, used by Bardeen, by his other 
gauge-invariant amplitude Vg, used by us. In our gauge, 
which is shear-free, dx'^/dt is equal to Bardeen's gauge- 
invariant Wg. 

The replacement of by Vg is necessary, because the 
field with the amplitude Vc, used by Bardeen, is directly 
related to the angular velocity of matter relative to spin 
axes of local gyroscopes at the sources throughout the 
universe. This is not a measurable input without prior 
knowledge of the solution of Einstein's equations, i.e. of 
goi, Ag, and Bg. — In contrast, the field with the am- 
plitude Vs, used by us, is directly related to the angular 
velocity of matter measured relative to geodesies (on Ef) 
which start at our position and go in the directions of the 
chosen spatial axes along our world line (resp geodesies 
to asymptotic quasars in an open asymptotically unper- 
turbed FRW universe). The input needed in the context 
of Mach's principle is the field with the amplitude Vs- — 
The difference between Bardeen's two gauge- invariant ve- 
locity amplitudes is 

(«c-«s) = -*. (17) 

The 5* term is a geometric output, hence it must be 
moved to the left-hand side, the geometric side (output 
side), of the momentum constraint. For the prefactors 



on Bardeen's right-hand side we use the "second" Fried- 
mann equation, 

(dH/dt) - ^ = -4^Gw(/5 + p). (18) 

Combining the first and second steps gives the new op- 
erator on the left-hand side (geometric output side) of 
the momentum constraint: [A -|- (2 + 2 — 4:)K/ — fi^] = 
[A — /z^], i.e. the three curvature terms cancel. The new 
term on the right-hand side (measured input side) is our 
energy current J^. This gives Eq. ([5|). 

From Eqs. dm [131) we obtain vl = dx'/dt = - A^,. 
Therefore the difference between the momentum den- 
sities (p + p)vl, used by Bardeen, and (p -I- p)vl, used 
by us, is analogous to the important difference between 
the canonical momentum pi = dL/dx^ and the kinetic 
momentum A:* = mdx''/dt in Lagrangian mechanics for 
point particles of charge q in an electromagnetic field, 

k' = p'-qA', (19) 
ip + p)vl - ip + p)[vl-Al], (20) 

where the 3-indices are raised and lowered by the 3- 
metric. The kinetic (resp. canonical) momentum is (resp. 
is not) measurable without prior knowledge of and Ag. 
The kinetic momentum density is equal to the LONB 
components T? = {p + p)dx~^/dt, while the canonical 
momentum is equal to the coordinate-basis components 

n = ip + p)glf [dx^/dt + A^]. 

To prove the postulate formulated by Mach (exact 
dragging of inertial axis directions) the kinetic momen- 
tum, resp the kinetic angular momentum of matter is 
the appropriate input. The canonical angular momentum 
(e.g. aroimd the z-axis: T^), cannot be used, because its 
measurement resp. definition needs an absolute element 
in the input. — More details on kinetic versus canonical 
momentum and angular momentum in Sec. VIII of [?[• 

A cosmological term in Einstein's equations is diagonal 
in a LONB, hence a cosmological constant cannot appear 
in the momentum constraint in LONB components. If a 
cosmological constant A is absorbed in p and p, one has 
(pa /pa) = — 1, and evidently A cannot make a contribu- 
tion to = (p + p)il. 



III. RIEMANN TENSOR UNIQUELY GIVEN 
BY OUR Bg FIELD 

In this section we show that our choice of FIDOs and 
hence our {Eg,Bg) fields are singled out uniquely. We 
specialize to gravitomagnetism on a Minkowski back- 
ground, since the material in this section will not be 
needed later, although it is very important: Our specific 
choice of FIDOs is singled out by the property that only 
with our choice of FIDOs the {Eg, i3g)-fields uniquely de- 
termine the Riemann tensor and the gravito-electric and 
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gravito-magnetic tidal tensors i3j-) for vector per- 
turbations, 

-Rm = (21) 

-^C^m^'lMRpm = ^^=Bl^y (22) 

Rki^^^rsRy.^, = ^(curl4)^. (23) 

The purely spatial Riemann tensor is zero for vector per- 
turbations, i.e. the intrinsic geometry of each slice St 
remains unperturbed, as explained at the beginning of 
Sec. lIBl — A vertical bar denotes the covariant derivative 
in 3-space, and a round bracket around two indices de- 
notes symmetrization. — The equation for is Eq. (35) 
of our companion paper 

The first equalities in Eqs. ([211 [22]) directly follow 
from the operational definition of the Riemann tensor: 
(1) i?j5jQ gives the relative acceleration of neighboring 
free- falling, quasistatic particles (geodesic deviation), i.e. 
the gravito- electric tidal field £i~-. The tidal 3-tensors 
{£lj,B~-j) must be traceless by definition. The 3-trace 
of RiQjQ is a 3-scalar, therefore it is zero in the vec- 
tor sector. — (2) The Hodge-star dual (— *^)io3o = 
-l^ippq^pm = ¥\pqRpm gives the relative precession 
of neighboring quasistatic gyroscopes' spins. The 3-trace 
of i^ipq-^pqjd) ^^^'^ because of the cycHc identity of the 
Riemann tensor, Rabcd + Racdb + Radbc = 0. Replacing the 
relative precession by the gravito-magnetic tidal field 
gives a factor {—\) according to Eq. @ for any choice of 
FIDOs. 

The second equalities in Eqs. (|2HI23p are only true for 
our (i?g, i3g)-fields, defined via our specific choice of FI- 
DOs. Our FIDOs are singled out by the property that all 
13 nonzero Riemann components are uniquely determined 
by spatial derivatives of our 2 vector fields {Eg,Bg). If 
one also uses time-derivatives of the Bg field and uses 
curling = —dtBg from Eq. (29) in [7], one concludes that 
our Bg field by itself determines all 13 nonzero compo- 
nents of the Riemann tensor: Our Bg field gives an ef- 
ficient and physically transparent way to represent the 
entire information in the Riemann tensor. 

The Weyl tensor Cap-yS (defined by having no non- 
trivial 4-contractions) and its Hodge dual *C, have been 
used in the literature |39l[40l.l4H to define the electric and 
magnetic tidal 3-tensors by contracting twice with a time- 
like unit vector. However, our tidal tensors have a differ- 
ent normalization, because tidal tensors must be directly 
connected to relative acceleration and relative precession 
(and hence to relative gravito-magnetic fields). The nor- 
malization of our tidal tensors in terms of the Weyl tensor 
is given by Qaja = and ^Qg^g = -i^-. The fac- 
tor i in the first equation comes from Cjg-g = 5^iojo 
valid for vector perturbations. The factor {—\) in the 
second equation comes from the conversion from the rel- 
ative precession (given by the Riemann tensor) to the 
relative gravito-magnetic field. 



The curvature invariant is R • R = RapjsR"^''^ = 
8[£ij£ij — jBijB~-j], and the curvature pseudoinvariant 

is *R • R = (iea/3p.i?%5)i?"^T* = - 4£.3.Ky.. These in- 
variants are not useful for first-order perturbation theory, 
because local Lorentz transformations of £:■-- and are 
second-order in perturbations. 

The connection in the opposite direction, finding our 
Bg field from the Ricci tensor: Any vector field in the 
3- vector sector, i.e. with zero divergence, is the sum of a 
field determined by its curl plus a harmonic field, which 
in Euclidean 3-space is a homogeneous field, because we 
require that the i?g-fields are not singular at spatial in- 
finity. Therefore the Ricci tensor determines our Bg-ficld 

up to a homogeneous Bg-field, which is equivalent to 
a time-dependent spatially rigid rotation of all FIDOs 
according to Eq. — Our solutions of the equations 
of cosmological gravitomagnetism on FRW backgrounds 
with K = (0, ±1), e.g. Eq. ([7]), are form-invariant un- 
der time-dependent spatially rigid rotations, see Sec. II HI 
The exact dragging of inertial axes, Eq. ([9]), is indepen- 
dent of the choice of FIDOs. 

Before Einstein's General Relativity, gravitomagnetic 
and gravitoelectric fields were discussed by Oliver Heav- 
iside '42'! based on a gravitational-electromagnetic anal- 
ogy. With his equations Heaviside could have derived 
the partial dragging of inertial axes by a rotating star 
(Lense-Thirring effect), although a factor 4 would have 
been missing. 



IV. EIGENFUNCTIONS OF THE LAPLACIAN 
FOR SCALAR FIELDS AND VORTICITY FIELDS 

IN (_B^s^i^^) 

Cosmological vorticity fields are of a different type 
from vorticity fields typical in fluid dynamics and con- 
densed matter physics with their vortex lines (line singu- 
larities) and fields often idealized with cylindrical symme- 
try. In cosmology we consider vorticity fields (= vector 
fields with vanishing divergence) which are regular, and 
we expand them in a spherical basis. These same fields 
also occur in the multipole expansion of electromagnetic 
radiation fields [i^. In this section we generalize these 
fields to the curved 3-spaces and H^. — In the first 
paragraph of Sec. lIVBj we explain, why the methods used 
in textbooks and in the literature for CMB-analyses are 
not adapted for our purpose. 

Our aim is to solve the momentum constraint with its 
operator (A — /i^), i.e. the inhomogeneous Helmholtz 
equation. As a preparation for the Green function of 
(A — we derive the eigenfunctions of the Laplacian 
for {E^, S^, H^) inside resp. outside a shell with sources. 
First we do this for scalar fields. 



10 



A. Eigenfunctions of the Laplacian for scalar fields 

in (£;^5^//^) 

Our new results in this subsection: We give the radial 
functions of the scalar harmonics in a simple, identical 
formioT all three geometries {E^, S^, H^), e.g. Eq. ([50]) . a 
form which is familiar from the spherical Bessel functions 
for elementary wave physics in Euclidean 3-space. This is 
in contrast to the complicated derivations and results in 
the literature. — In Subsec. flV Dl we shall show that the 
radial functions are also identical for scalar and toroidal 
vorticity harmonics, if and only if one uses those vector 
spherical harmonics, X^^^^ which have point- wise norm 
and LONB components independent of R for fixed {9, (f). 

Since the momentum constraint Eq. ([5]) is an ellip- 
tic equation (no partial time-derivatives), we consider, 
in this and all following sections, a FRW universe with 
ii' = ±1 at a given time with units of length chosen such 
that the curvature radius is Cc = a = 1 at the given time. 
We also include the case of FRW with K = (Euclidean 
3-space E'^) at a given time with the scale factor a = 1, 

a = 1 (24) 
ds^ = dx^ +RK{x)[dO'^ +s^^^Od(j)'^] (25) 
Rk{x) = {X, sinx, sinhx}, K = {0, +1, -1}. (26) 

The following relations will be useful, 

R" = -KR, R'^ = 1- KR^, (27) 

where ' = d/dx, and R is short for Rk{x)- 

For scalar fields the angular eigenfunctions of the 
Laplacian A are spherical harmonics Yfm(0, (/>) with 
eigenvalues —£{£ + 1)/R^. — The radial eigenfunctions 
of the Laplacian are the generalization to A' = ±1 of 
the spherical Bessel functions jiiqr), which are the ra- 
dial eigenfunctions in Euclidean 3-space. These general- 
ized spherical Bessel functions, the hyperspherical Bessel 
functions, will be denoted in this paper by j^^^ (x) to em- 
phasize the close correspondence to Jgf~°\x) = je{qx)- 
We shall always denote the eigenvalues of the de Rham- 
Hodge Laplacian A by (— fc^). For A acting on scalar 
eigenfunctions and on vector or tensor eigenfunctions in 
the scalar sector we define q by 

scalar sector: q'^ = + K. (28) 

It will turn out that q is the wave number in the radial 
oscillations sin(gx) in Eqs. dSOl ED EH) • 

The radial part of the Laplace operator acting on a 
scalar function f{x) is A/ = R-^ d_^R'^ d_jy y^j^g 

R-^-^R = ~K we obtain A/ = R-^^{Rf) + Kf. 
Hence the radial eigenfunctions for scalar harmonics sat- 
isfy 

Rd^(^3,i )-[-<i ■ (29) 



The solutions regular at the origin are 

which we call generalized spherical Bessel functions. 
Writing the radial equation and the solutions in this 
form gives a unified notation for all three geometries 
{E^, S^, H^). Compared to spherical Bessel functions (fa- 
miliar from wave physics in Euclidean 3-space) the only 
changes are the replacements of the independent variable, 
r — > X, in the radial oscillations and in the derivative, 
and the replacement of the denominators and pref actors, 
r — > R{x)- — From Eq. ([50)1 we obtain the generalized 
spherical Bessel functions fq^\x) for I — (0, 1), 

jljlloix) = ^sin(gx), (31) 

iiEiW = -J2^^^^^^ix)]- (32) 

The generalized spherical Neumann functions are ob- 
tained from Eqs. ([501 EU EH) by the replacements 
sin((7x) ~ cos(gx) and cos(qx) ~^ +sin((7x). The gen- 
eralized spherical Hankel functions are given by /i^^'^-* = 
j ± in. — Relevant for Mach's principle are the spheri- 
cal Bessel and Hankel functions for I — 1 generalized to 
K ^±1. They are given by Eq. (EH) and 

KSf^^x) ^±^^^^[]^ exp(±zgx)]. (33) 

The eigenfunctions of the Laplacian will be used in 
the context of Green functions with sources at x = Xsi 
hence they will be used only for < x < Xs resp. for 
X > Xs- Therefore there are no geometric restrictions on 
the eigenvalues in our context. 

The simplicity of the elementary functions p2i I33p . 
their identical form for all three geometries, the familiar 
appearance from elementary wave physics in Euclidean 
3-space, and the straightforward derivation are in con- 
trast to the complicated derivations for the hyperspher- 
ical Bessel functions $£(x) via the associated Legendre 
functions v) discussed extensively in [i^l and 

used by [11] . 

B. Vector spherical harmonics 

In this subsection we introduce the 3-dimensional vec- 
tor spherical harmonics {X^^, CyYirn) and the Regge- 
Wheeler harmonics {x^^,d^Yern), and we obtain their 
curls and divergences in Eqs. (jH]) - ([^5]) . which are valid 
and form- identical not only for the 3-spaces {S^,H^, E^), 
but for any spherically symmetric Riemannian 3-space. 
Our curl- and div-equations for this general case are very 
much simpler than the equations in the literature for the 
special case of Euclidean 3-space. — The general and 
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simple equations (|44p - ([49]) are all that 's needed to con- 
struct vorticity fields (divergenceless) and to explicitly 
write down the de Rham-Hodge Laplacian on vorticity 
fields, AV — — curl curl V, for any Riemannian 3-space 
with spherical symmetry. — Textbooks usually present 
the vector spherical harmonics of [i^, 113], which are 
not adapted to cosmological perturbation theory, because 
they mix the vector and the scalar sector, see the end of 
this subsection. — The spin-s spherical harmonics used 
in CMB analyses ^l| are not adapted to our problem, 
because they mix parities, i.e. the toroidal and poloidal 
sectors of vector spherical harmonics, see the end of this 
subsection and Sec. IIVCI 

Vector spherical harmonics X^^{0, 4>) form a basis for 
vector fields tangent to the 2-sphere. In the notation 
of differential forms, the vector spherical harmonics of 
Regge and Wheeler [48|, which we denote by a lower case 



are given by 



^tm — 



dYi, 



im ; 



(34) 



i.e. the gradient of Y^m, and its Hodge dual on the 2- 
sphere, * . In component notation Eq. ([M)) gives 



da Ye 

m 1 



1/3 g'^^ d^Ye^, 



where a — 
combination e 



(35) 
(36) 

EE +1. In the 



fj,i,g'^ = the radius R of the 2-sphere 

drops out, Eg"^ = (sin6')^-^, e^^ — — (sin0). 

The vector harmonics of Regge and Wheeler have co- 
variant components (1-form components) which are in- 
dependent of the radial coordinate Xj 



(37) 



In contrast, the 'physical' vector spherical harmonics 
are denoted by the upper case X^. They are used in 
the literature on classical electrodynamics [4^ and called 
pure-spin vector harmonics by Thorne [ioj . They are de- 
fined with an i?-factor which makes the radial covariant 
derivative vanish. 



0. 



(38) 
(39) 



This is equivalent to requiring that the point-wise norm 

di-^im' ^ini) ^^'^ LONB components of are in- 
dependent of R for fixed (0, 0). 

We rewrite the 2-dimensional Hodge dual of Eq. 
as * VYfm = — x VYgm- The generator of rotations, 
the angular momentum operator of wave mechanics in 
units of ?i, is L = —i{R x V), where R = R{x)e^. Hence 



Ylm, 



RVYtrr 

iLYfm 



(40) 



{R X v)r,, 



?X X Kn (41) 



(42) 



where we have also listed the third basis field needed for 
3 dimensions, e^Yim- The triple (X+, X~, e^) is point- 
wise orthogonal and has positive orientation. X+ and 
X~ have the same norm point-wise. 

The parity is P — (^1)^ for X^^^ and e-^Yim: while 
P= (-1)^+1 for 

All three basis fields {X^,e^Ygm} are eigenfunctions 
of the total angular momentum operators {J^,Jz} with 
eigenvalues {({i -\- l),m}. This follows because exterior 
differentiation d, taking the Hodge dual (*), and mul- 
tiplication with commute with rotations of the total 
system. 

The normalization and orthogonality relation on the 
2-sphere at any R is 



I'm'' ^ I ^^^Im^^l'm' ^ 

= i{t + l)6ut5mm'Spp' ■ (43) 



This relation is crucial for the projection of a general 
velocity field on the term with —l^P — -\-). The global 
(Hilbert space) scalar product on the 2-sphere is denoted 
by (y,W), and <V,W> denotes the point-wise scalar 
product g{V* , W), where the superscript * stands for the 
complex conjugate. 

Computations of curl and div and results are much sim- 
pler for Regge- Wheeler harmonics than for the 'phys- 
ical' harmonics X^. — Also the basis fields [e^R~'^Yem] 
are computationally much simpler than the basis fields 
[e^Y^m], because the divergence of the former vanishes 
apart from a (^-function at the origin, which turns out to 
be irrelevant in our context. 



diva; 



curl: 
div; 



'tm 



1(1 
i?2 



Yt 



curl X 



i?2 



div [Yi„^^] 

curl [Yem^] 



= 



i?2- 



(44) 

(45) 
(46) 

(47) 
(48) 
(49) 



If we had computed the divergence and curl of Xf^ and 
Ygmex^ there would have been twice as many terms, and 
the terms would have been more complicated. 

Computations of curlV^ = WxV and divV = W-V, 
where V is the covariant derivative, are made very simple 
with the calculus of differential forms, where no Christof- 
fel symbols are needed, see Eqs. (|A.21 IA.5]) . 

The spin-s spherical harmonics of Newman and Pen- 
rose 0, HO, ISH, which are widely used for cosmic mi- 
crowave anisotropics f3l'|, are not adapted to our prob- 
lem, because they mix parities, i.e. they mix toroidal and 
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poloidal vorticity fields (see Sec. IIV Cp . while the drag- 
ging of incrtial axis-directions is caused by energy-flows 
with 1+, i.e. in the toroidal sector (see Sec. IV A|l . 

In most textbooks vector spherical harmonics are con- 
structed by coupling a definite orbital angular momentum 
with spin-1 basis vectors to obtain a definite total ang ular 
momentum via Clebsch-Gordan coefficients l46l 14711491. 
These basis states are not adapted to cosmological per- 
turbation theory, because they mix the poloidal vorticity 
sector and the scalar sector^ see the following subsection. 



C. Toroidal and poloidal vorticity fields 



For 3-dimensional vorticity fields the 3-divergence 
must be zero. The vector spherical harmonics for 
TO = are vector fields purely along the direction of 
± 60 and independent of cf). Therefore the x'[^ are named 
'toroidal and evidently their divergence is automatically 
zero. The toroidal harmonics with m ^ Q are generated 
from those with m = by rotations, and their divergence 
is again zero, as also seen from Eq. (|46p . — Multiplica- 
tion with a radial function 17*°' (x) gives the 3-dimensional 
toroidal vorticity fields. 



(50) 



The 3-divergence is not zero for the vector spherical 
harmonics x'^^ by themselves, divx^'J^ 7^ 0. To obtain 
vector fields V with zero 3-divergence we must add a 
part perpendicular to 5^, e'^F{x,0,4>), where F can 
be expanded in scalar spherical harmonics Y^m. Hence 
the 3-dimensional poloidal vorticity fields can be written 



VZ\x.0.c^)=gf°''''^x 



(51) 

Both terms have P = (—1)^. For to = this 
equation gives vorticity fields V with only radial and 
0— components, since x^^^^ points along the meridians. 
Hence the name 'poloidal' for these vorticity fields. 

The condition for vanishing divergence of the poloidal 
vorticity fields is evaluated using the vector identity 
div (gV) = g div V + V ■ grad^ and Eqs. 



(pol,rad) 



dx 



£(^+1)5) 



(pol.tg) 



(52) 



Again, the equations in this subsection are valid for any 
Riemannian 3-space with spherical symmetry. 



D. The radial functions for vorticity harmonics 

Relevant for rotational frame-dragging are energy cur- 
rents and vector potentials in the toroidal sector as 
demonstrated in Sec. IV Al Hence we now focus on the 
toroidal eigenfunctions of the Laplacian. — Poloidal en- 
ergy currents cannot contribute to the precession of a 



gyroscope at the origin. Therefore poloidal eigenfunc- 
tions of the Laplacian are not considered in this paper, 
but they will be needed and presented in a forthcoming 
paper, "The other half of Mach's principle: Frame drag- 
ging for linear acceleration" ,52,]. 

In this subsection we derive the following important 
and simple result: 

1. In the eigenfunctions of the Laplacian for toroidal 
vorticity fields on {S^ , , E^): The physical vector 
spherical harmonics are multiplied by radial 
functions which are identical to those for scalar har- 
monics, i.e. the generalized spherical Bessel func- 
tions 2]^^ (x) of Eq. ([50]) resp. the corresponding 
Neumann and Hankcl functions. 

2. For toroidal vorticity harmonics we have the rela- 
tion q = k: 

(A-fg^)[4f (X)X,„(0,0)]=O. (53) 

In all sectors q denotes the wave number in the 
radial oscillations sm{qx), and (— fc^) denotes the 
eigenvalue of the de Rham-Hodge Laplacian. Note 
that the relation for scalar harmonics is q^ — 
P+K, Eq. (ESI). 

To prove this result, all that's needed is the expres- 
sion for the de Rham-Hodge Laplacian for vorticity fields, 
A = -curl curM, and the simple Eqs. dH]) - (gH]). — For 
the rest of this subsection we drop the subscripts {i, to). 
We start with a toroidal vorticity field A*™ = g^"^^ ^-^d 
evaluate its curl. Using the vector identity curl (gif) — 
5 curl u -I- (gradg) x v and Eq. (|17|) gives 

curl A*" = cm\{g'^'x~) 

= -i-^aT)^^ - m + i)gT[YR-^e^]- (54) 

In a second step we take the curl a second time, which 
gives minus the de Rham-Hodge Laplacian of 



-curl curl A 



i?2 9a\x 



(55) 



where we have used Eqs. (|45l [49]) . Up to here all for- 
mulae in Sees. llVB] - IIV Dl are valid for any Riemannian 
3-space with spherical symmetry. We now specialize to 
{S^ , , E^). We rewrite the toroidal vorticity fields in 
terms of the physical vector harmonics X^^^ = Rx^ 
and their associated radial functions G'*°'' = R~^g*^°'^, 
hence g^°'^x~ = G*^™X~. Comparing Eq. (|55p rewritten 
in terms of G'°'' with Eq. ([29l) gives the main result of 
this section, Eq. (l53|) . 

In contrast to ^16] we have given a simple, unified treat- 
ment of all three cases, K = (0, ±1) , with form-identical 
derivations and results. — Ref. [l6| discussed the solu- 
tions of the homogeneous momentum constraint (which 
are eigenfunctions of the Laplacian) without referring to 
the Laplacian. 
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V. MACH'S PRINCIPLE 

A. Symmetries and Mach's principle 

We treat all linear perturbation fields on all types 
of FRW backgrounds, K — (±1,0), and all energy- 
momentum-stress tensors, i.e. all types of matter (also 
dark energy and a cosmological constant), not necessarily 
of the perfect-fluid form, and for totally general field con- 
figurations of energy currents J| = T'i . This is in striking 
contrast to the previous literature, which only treated the 
artificial situation of spherical shells of matter rotating 
rigidly around one given axis. — However the mathemat- 
ics of totally general perturbations can be reduced to the 
mathematics of the special case of spherical shells of mat- 
ter in rigid rotation around a given gyroscope with the 
help of two theorems, which are based on the symmetries 
relevant for Mach's principle, rotation and parity: 

1. The precession of a gyroscope (relative to given local 
LONB axes) cannot be caused by scalar perturba- 
tions nor, in linear perturbation theory, by tensor 
perturbations. — In the vector sector the preces- 
sion of a gyroscope can be caused only by energy- 
current fields with — 1"*" relative to the given 
gyroscope's position, i.e. of the toroidal type and 
with i=l. 

2. On every mathematical spherical shell centered on 
the gyroscope considered: The field component (in 
the harmonic decomposition of the general energy- 
current field) which is toroidal and has £ — 1 (rel- 
ative to the gyroscope) is given by an equivalent 
rigid rotation with an equivalent angular velocity of 
matter 5l™^"^°''(xs)- The equivalent angular veloc- 
ity of matter is given by the global scalar product of 
the Jg-field with the toroidal fields X^^^ ^ on the 
shell with radius x ■ 

Jdn< XiJ, „^ , Ux, 0, 0) >^ (Je)7=i.™(x) 

= - V/1W3(P0 +P0) R{X) [^ra{x)Y::l'\ (56) 

where < ....... > denotes the point- wise inner 

product, and dfi is the element of solid angle, while 
f2m denotes spherical-basis components of the an- 
gular velocity. The transformation between the 
spherical and the Cartesian components of Vl is 
given by Vlm=o = riz, and ^m=±i = {t^x - 
iVly)/y/2. 

The proof of theorems (1) uses the following facts: The 
precession of a gyroscope relative to given LONB axes, 
d{Si)/dt, which equals the torque on the gyroscope, is an 
axial vector, = 1"*". The precession of a gyroscope is 
caused locally by the gravitomagnetic field -Bg(O), which 
is also an axial vector. — For scalar perturbations all 
fields are derived from scalar fields (via differentiation). 



but this can only produce source-fields in the natural 
parity sequence, 0"*", 1~, 2"*", etc, which cannot contribute 
to the precession. — For tensor perturbations (gravita- 
tional waves) all linear perturbations are given by a trace- 
less, divergenceless 3-tensor, from which one cannot form 
an axial vector field to generate a Bg field at the origin. 

The proof of theorem (2) uses the following facts: For a 
general energy-momentum-stress tensor (not necessarily 
of the perfect-fluid type) the local center-of-mass veloc- 
ity V is given by = (po ~^ Po)~^'^1 iii linear pertur- 
bation theory, where po and po refer to the unperturbed 
FRW background. — For fixed X: ^ toroidal velocity field 
^fci m is a fiow equivalent to a rigid rotation of a shell 
of matter. 

The gravitomagnetic vector potentials Ag relevant for 
the gyroscope's precession, must also be toroidal with 
£ — 1. The gravitomagnetic field Bg — curMg, which 
causes the gyroscope's precession, must be poloidal with 

B. Green functions for (A — /i^) on vorticity fields 
in open FRW universes 

We derive the Green function for the gravitational vec- 
tor potential Ag for K = (—1,0) generated by an energy- 
current source which is toroidal, has £ = 1, m = 0, and is 
concentrated at the geodesic distance Xs relative to the 
gyroscope considered. Such a source corresponds to a 
shell of cosmological matter in rigid rotation around the 
z-axis, 

v'^^n,S{x-Xs). (57) 

The Green function for arbitrary £ is not needed to 
compute the dragging of inertial axes at the origin. The 
derivation for general £ is totally analogous to the case 
£ = 1 treated here, and it gives the general solution of 
the momentum constraint for the toroidal sector. 

As noted at the beginning of Sec. IIV A[ all equations 
in Sees. HVl and IVl involve the perturbation quantities at 
one time, and at the chosen time we set the scale factor 
a equal to one. 

Inside and outside the rotating shell, A is a toroidal 
vector eigenfield of (A — /i^) according to the momentum 
constraint Eq. ([5]). The wave number q in the radial 
oscillations is g = ±i/i according to Eq. ([55)1 . We choose 
the positive sign. 

Inside the shell the radial function G{x) — fix) which 
multiplies X^^-^ according to Eq. is the generalized 
spherical Bessel function for K — (±1,0), Eq. (|32p . 

»A^'j^i,,=i(x)-^[|sinh(^x)]. (58) 

Note that any regular toroidal vector field must vanish at 
the origin. 

Outside the shell the radial function for an open uni- 
verse, K = (—1,0), which is regular at spatial infinity 
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for q = ifjL, is the generalized Hankel function of the first 
kind, Eq. (|33p . exponentially decreasing for x — > cx), 



~d 1 

'dxR 



exp(-^x)] (59) 
— Yukawa force. (60) 



Note that the Yukawa potential is (1/i?) exp(— /ix), while 
the Yukawa force is minus the gradient (the x-derivative) 
of the potential. — For a closed universe this solution is 
not acceptable, because the 1 /R factor makes it singular 
at X = Ti", the antipodal point to the origin. This will be 
treated in subsection IV Dj ^ 

The vector potential is tangential to the shell, 

and it must be continuous across the shell, because the 
component of B normal to the shell must be continuous, 
otherwise div_B 7^ on the shell. In an open universe, 
K = (—1,0), we obtain 



fT{x,Xs) = Nj{x<)h^'\x>), 



(61) 



where we have dropped the subscripts {q = ifi, £ = 1} 

and the superscript {K — —1,0) for better readability. 

The arguments x< resp. x> denote the smaller resp. 

larger of the two arguments (XiXs), and Xs refers to the 

radius of the shell. 

The normalization factor is determined by the 

strength of the source Je on the shell: Einstein's G"- 

4, 

equation, Eq. ([5]), determines the discontinuity of the 
tangential component Bg, 

(curl B = Six - Xs) disc[S,-] = -16^Gn J| 

= -(5(x-Xs)16^GN(p + p)ai?sSin0. (62) 

The discontinuity of Bg is determined by the discontinu- 
ity of its radial function For an open universe, 
K^i~l,0), Eqs. dSlElinil) give 

disc [/(P°'^*^)] = -disc[A/^o..] 

= -Ndisc[-f{j{x<)h(^\x>)}]s 
dx 

= -NW[j,h('^U, 



(63) 



where we denote the Wronskian by W[f, g] = [fg' — f'g) 
with / = Our Wronskian multiplied with is inde- 
pendent of X, as can be seen from Eq. Therefore we 
evaluate R'^W in the limit x ~^ 0. We use j^^l^ e=i(x) ~^ 

{^X^^/3){l+K/fi^) andn(f]^^,^i(x) ^ (/xx)"" forx « 1, 
and we obtain 



I K 

(1 + -) (64) 



valid for K — (±1,0) and for all R. We note that for 
l,m = 0) 



and all other components are zero. We multiply 
'-i=i m=o)e to obtain disc[i?g], which 



disc[/]P°''*^'] with iX+ 



we insert in Eq. (p^ to determine N, 



K 



N = -{\&^G^){p + p)Rln,^i{\ + — )-iy4^. (66) 

This completes the computation of the Green function, 
Eq. (|6ip . in an open universe, K = (—1,0), for the vec- 
tor potential generated by the toroidal current with 
£ = 1, m = at Xs, specifically by the velocity field of 

Eq. dnzi). 

The gravitomagnetic field B at the origin, the posi- 
tion of the gyroscope, is obtained from Eq. Bg{0) 
written for arbitrary orientation of fig is 



BiO) = -217gy„(0) - -ns^[167TGNip + p)]R'sY^ixs), 

(67) 

valid in an open universe, K ~ (—1,0). This equation 
shows, how (is, the angular velocity of the rotating shell, 
determines r^gyro, the angular velocity of precession of a 
gyroscope at the center of the rotating shell. 
In Eq. ((67)) we denote the prefactor by fP, 



~2 



167rGN(p + p) 



2 = -4{dH/dt). 



(68) 
(69) 



The connection between fl and /i is given by the (dH/dt) 
equation, which follows by subtracting the first from the 
second Friedmann equation. 



(dH/dt) 



= 1 => 



^2' 



K 



K 



= -47rGN(p+p) 
- (-)'• 



(70) 



(^ti,™=o)<^ = (^7=i,™=o)^ = -Vm^)^^^0, (65) 



Three length scales occur in our problem: the -ff-dot ra- 
dius (^/2)~^, the (p -|- p)-radius (/i/2)~^, and the curva- 
ture radius Cc = 1. 



C. Exact dragging of inertial axes for 
perturbations of an open FRW universe 

The main result of this paper: (1) The equations for 
"Dragging of Inertial Axes", Eqs. ^ [TU|) . are obtained 
directly from Eq. (j67p by integrating over all shell radii. 
(2) The equation, which shows that there is exact drag- 
ging (as opposed to partial dragging) of local inertial 
axis directions, is the statement that the weight function 
W^(x) has its normalization equal to unity, Eq. (|lip . 

We now give the proof that the normalization of the 
weight function is indeed equal to 1 for the case of an 
open universe, K = (—1,0). First one does a partial inte- 
gration of [d/dx) in the Yukawa force. This gives Wtot — 
p,^ J dxRR' e:xp{—ij,x)- A partial integration of the factor 
exp(-^x) gives Wtot = / dx{l~2KR^)pL-^ exp(-/zx), 
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where we have used Eq. ((27l) . If instead of the last step, 
one performs a partial integration of the factor {RR'), 

one obtains Wtot = P-^ J dx{R^ / '2') IJ. exp{— ^x) ■ Adding 
/i^ times the first expression plus {4:K) times the second 
expression one obtains {iK + /i^)Wtot — A^, and using 
Eq. m gives w/j^^"'^"^ = 1. 

Our answer to Mach's question "what share?" is given 
by (1) the radial weight function W{x) with its Yukawa- 
force cutoff, and (2) the angular weight function X^J[ ^ 
of Eq. in the projection of the general velocity field 
on the toroidal vorticity sector with £ = 1. 

The kinetic angular momentum L per dx is directly 
obtained from rJmatter using the moment of inertia per 
dx of a fluid shell, which is (87r/3)(p + p)R'^, 

^ = Y(P+P)i?'f^matter(x). (71) 

Evidently the kinetic angular momentum L is determined 
by angular velocity measurements without prior knowl- 
edge of 500; which is the output of solving Einstein's equa- 
tions. 

The gravito-magnetic moment per dx, dflg/dx, is de- 
fined in the same way as in ordinary magnetism, except 
that the current of charge is replaced by the current of 
energy X : 

j dn[nx J,]ra = -Y J dn< X^J^^^ ■ J, >,(73) 

where n is the radial unit vector at the source point, and 
R = Rn. The last two equations show that the gravit- 
omagnetic moment and the angular momentum involve 
a projection of the energy-current field on the sector of 
toroidal vorticity fields with i = 1. 

rjgyro is determined by an integral over the density of 
kinetic angular momentum: From Eqs. ((9l I7ip 

4(P) = -2n,y,, = -4Gn J d(volQ) 

(74) 

Expressing the kinetic angular momentum L by the 
energy current we obtain the fundamental law for cos- 
mological gravitomagnetism in integral form for an open 
FRW universe, (if = -1, 0) : 

Bg(P) = -2dgy,.o(P) = 

= -4Gn j d{vo\Q)[npQ X ■UQWf.iXPQ)- (75) 

Again: figyro is given by the sources Jg at the same time. 
This follows directly from: (1) The facts that dragging 
effects are given directly by the momentum constraint, 
and that the momentum constraint in the vorticity sec- 
tor, Eq. ([5]), for the general time- dependent context of 



gravitomagnetodynamics is an elliptic equation (it has 
no partial time-derivatives). These facts have already 
been recognized by [lll,[l3,l3| among others. — (2) The 
fact that the tensor sector (gravitational waves) in lin- 
ear perturbation theory (and of course the scalar sector) 
cannot produce a gravitomagnetic field i?g nor a torque 
on a gyroscope has been demonstrated in Sec. IV Al 

The vector structure oiJLqs. (|721I731 [75 | needs explana- 
tion, since there is no global parallelism on a hyperbolic 
3-space . This is explained at the end of Sec. II Gl 

The differences between Eq. ()75|) and Ampere's law in 
integral form are: 

1. the replacement of the current of electric charge Jq 
by the energy current J^, 

2. the factor (— Gn), as in the transition from 
Coulomb's law to Newton's law, 

3. the additional factor 4, which occurs in the transi- 
tion from Ampere's law of ordinary magnetism to 
gravitomagnetism, 

4. the replacement of the 1/r^ force in Ampere's law 
by the Yukawa force (%) for cosmological gravit- 
omagnetism, 

5. the need to distinguish between R (in the denomi- 
nator of the Yukawa force) and x (in ttie exponent 
and in the derivative) for FRW backgrounds with 
curved 3-space, K = ±1, 

6. the remarkable fact that the law (|75|) for cosmo- 
logical gravitomagnetism is valid for general non- 
stationary situations, while Ampere's law is valid 
only for stationary currents. 

What distance interval x of sources dominates drag- 
ging of local inertial frames at P? — For a spatially flat 
universe, the weight function W{x) grows linearly from 
X = to X « = A^^i and it decays exponentially 
for X ^ A*"^- Hence the sources Q around xpq = 
dominate the frame-dragging at P. — For a spatially hy- 
perbolic universe, K = —1, and for // 3> = 1 the 
situation is approximately the same as for K = 0, since 
the frame dragging is caused by sources at distances much 
smaller than the curvature radius. — But for K = —1 
and for x much larger than the curvature radius Oc = 1, 
the effective exponential cutoff factor in the weight func- 
tion W{x) is exp[— (^ — 2)x]. From {p-\-p) > and from 
Eqs. (El EOl) with K = -1, it follows that > and 
^ > 2. For fl « 1 we have (^ — 2) w /i'^/4, the cutoff 
factor becomes exp(— /i^x/4), and we have three regimes 
for W{x) ■ For < x << 1 the weight function W{x) 
grows linearly, for 1 << x << the weight function is 
constant, and for x >> it decays exponentially. On 
a logarithmic scale the sources near x = dominate 
frame dragging for /i << 1. — For all open universes 
the dragging of inertial frames is dominated by sources 
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at a distance equal or larger than the (p + p)-radius: Lo- 
cal inertial frames are rigidly "in the grip of the distant 
universe" . 

Form-invariance of the solution Eq. (|75p under trans- 
formation to a rigidly rotating frame: If the reference 
FIDO at the position of the gyroscope changes his lo- 
cal spatial axes to new ones, e?, with angular velocity 

— il*(0) relative to the old ones, Eq. (|75p remains valid 
as it stands, because both sides of Eq. ([75]) change by the 
same term. On the right-hand-side (n x J^) changes by 
[— [p+p) ri* (0) i? sin 9 eg] , where we have put the z-axis of 

the spherical coordinates in the direction of — The in- 
tegration over the solid angle gives [S-k /'i){p+p)R^*{Q). 
The integration over x is the same as in Eq. (dJ). The re- 
sulting change on the right-hand side is [— 251*(0)], which 
is equal to the change on the left-hand side. 

To obtain definite values for measurements of fields 
and i3g, one must make an arbitrary choice of a definite 
state of rotation of FIDO-axis directions along one world 
line, as explained at the end of Sec. II Bl But for any such 
choice: The field of measured energy currents all by 
itself completely determines the gravitomagnetic field i?g 
at all points, i.e. the precession ilgyro of gyroscopes at 
all points: No absolute element is needed in the input, as 
required in the starting point of Mach 0, Q . 

For cosmological gravitomagnetism and for given 
sources J^, measured via angular velocities, distances, 
and {p + p): The solution Eq. ([75]) for the momentum 
constraint is unique. There exists no ambiguity of adding 
a solution of the homogeneous equation. No regular so- 
lution exists for the homogeneous momentum constraint 
in cosmological gravitomagnetism, (A — /i^)Ag = 0. — 
Therefore for regular solutions no boundary conditions 
are needed at spatial infinity for open universes. See also 
Sect. XI of 0. 



D. Exact dragging of inertial axes for 
perturbations of a closed FRW universe 

To find the Green function for in a closed universe 
and in the toroidal sector, we note that a toroidal vector 
field V(Li TO which is a regular vector field at % = and 
at X = must have a point-wise norm which vanishes 
linearly at these two points. Therefore the radial solution 
= Pa^^ which multiplies X^^^ must also vanish 
linearly at x = and at x = tt. Hence the radial part of 
the Green function can be written as 



NW [ j{x),ji^-x)\ 



(77) 



fT{x.Xs)=Njix<)j{7r-X>)- 



(76) 



We have again dropped the subscripts {q — ip, £ — 1) 
and the superscript (if = +1) for better readability. 

The discontinuity in the derivative of the radial part 
of the Green function, the analog of Eq. ([55)1 . is 



disc 



d 



rtor ■ 
JA . 



To evaluate the Wronskian W we use the method given 
before Eq. ([M)) . Specifically, (R^W) is independent of 
X, and we evaluate this expression for x ^ 0- We use 

jSSiW -> (*/3)xm(1 + 1/m') ^ndjl^g7M^-x) - 
i/ {x^p^) sinh(/x7r) for x ^ 1, and we obtain 



w[j',V^x):j\y{^~x)] 



3W 



I K=+l 
q—ifi^ i—1 



Hence the normalization N is 



N = Nsinh^^ipn), 



(78) 



(79) 



where N = iv(^-" -i) is given by Eq. ((66)) . 

The gravitomagnetic field B at the origin generated 
by a rotating shell ri(x) = ^(x ~ Xs)f^s on a FRW back- 
ground with _ftr = -1-1 is given by taking the result ([67]) 
for K — (0, —1) and making the replacement 



exp(— /ix) =^ sinh ^{pTr) sinh[/^(7r — x)]- 



(80) 



For an arbitrary velocity field it directly follows with this 
replacement from the result of Eq. ([9]) that ^gym is equal 
to the weighted average of fimattcr , 



fJgyi-O 



W{x) 



dx f^ii 



r(x) W{X), 



-p R-^smh ipn) 
3 

^{isinhK--x)]}. 



(81) 



Again, the averaging weight function is normalized to 
unity, Wtot — Jq dxW{x) — 1- The proof uses the 
same steps as the ones given to prove unit normaliza- 
tion of W^(x) for open universes in the second paragraph 
of Sec. IV CI Hence for perturbations of a closed FRW 
universe the evolution of inertial axes exactly follows the 
weighted average of cosmic matter motion, there is exact 
dragging of inertial axes by cosmic energy currents, as 
stated in Mach's postulate. 

The fundamental law for gravitomagnetism in inte- 
grated form follows with the replacement (|80p from 
Eq. dZSl), 



BiP) 



4Gjvsinh ^/itt) / d{vo\Q) {npQ x X 



-^{-^[sinh[/i(7r-x)]}. 



(82) 



The solutions ([5T1 [5^ are again frame-invariant, i.e. 
form-invariant when going to a reference frame which is 
in globally rigid rotation relative to the previous reference 
frame. The proof is identical to the one given for open 
universes in Sec. IV CI 
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Einstein, Wheeler, and others, recently e.g. Bicak et 
al jl^], have held the view that from the Machian point 
of view closed universes are preferable. In contrast our 
results show that Mach's hypothesis, exact dragging of 
inertial axes, holds for linear perturbations of all FRW 
universes, open or closed. 

E. Exact dragging of inertial axes for 
perturbations of Einstein's static, closed universe 

For a static, closed universe /i^ = —4:{dH/dt) = 0, and 
(/i/2)2 ^ K ^ +1. Taking the hmit ^ ^ of the result 
for a FRW universe with K = +1, Eq. (15^ . one obtains 
the fundamental law for gravitomagnetism for Einstein's 
static, closed universe, 

/d 1 
d(volQ) (npQ X J,) - XA)](83) 

Again rigyro is equal to the weighted average of rimattcr, 

f^gyro ^ dx f^mattor(x) W (x) , 

Jo 

W{x) = (84) 

where Wtot = Jq dxW{x) = 1- Eq. states that also 
for perturbations of Einstein's static closed FRW uni- 
verse the evolution of inertial axes exactly follows the 
weighted average of cosmic matter motion, there is exact 
dragging of inertial axes by cosmic energy currents, as 
stated in Mach's postulate. — The solution ([55)1 is again 
form-invariant when going to a reference frame which is 
in globally rigid rotation relative to the previous refer- 
ence frame. 

Einstein's static, closed universe (3-spherc) has been 
important in the work of Ozsvath and Schiicking [s^, 
who added a Bianchi IX amplitude (the lowest mode of 
tensor perturbations) to Einstein's universe, and then in- 
vestigated Mach's principle. Unfortunately Ozsvath and 
Schiicking considered the vanishing of the local vorticity 
relative to a gyroscope to be a test of Mach's princi- 
ple. This is in contradiction to our conclusion that the 
weighted cosmic average of the energy currents in Eq. (j84p 
directly determines the time-evolution of gyroscope axes. 
Our result agrees with Mach's hypothesis that some (to 
him unknown) cosmic average of matter motion would 
determine the time-evolution of gyroscope axes. 

F. Mach's principle for the limits to the Milne and 
the de Sitter universe 

Many physicists (but definitely not Mach) have asked 
again and again: "If you could remove a body far away 
from all sources of gravity, what would be the inertial 



behaviour of this body?" — If the density of galaxies 
was e.g. 10~^ per Hubble volume in an otherwise empty 
universe, a particle could be a distance of 10"^ Hubble 
radii away from all other matter. All the same, if galax- 
ies in the unperturbed universe were homogeneously and 
isotropically distributed, we have shown that for a uni- 
verse with linear vorticity perturbations Mach 's principle 
still holds. 

The Milne universe is the limit p/pcrit ~* (for p/p 
fixed) of a FRW universe. The first Friedmann equa- 
tion gives H"^ + K/al 0, hence we have a hyperbolic 
universe with a curvature radius ac equal to the Hubble 
radius H~^. The Milne universe is equivalent to a part of 
Minkowski space-time, the forward light-cone originating 
at the space-time point where the Hubble expansion (of 
test particles at rest in the Milne universe) started. 

How can exact dragging of inertial axes by cosmic en- 
ergy flows work for p/pcrit ^ 1, when the total energy 
within a Hubble volume is arbitrarily small? According 
to the discussion in Sec. IV Cl the weight function becomes 
W{x) ~ (mV4) exp(-/i2x/4) for fl^ = IGnGt^ip + p) « 
1. The prefactor [fi? /A) is arbitrarily close to zero, but the 
integral for Wtot would be linearly divergent without the 
cutoff at Xcutoff — Afl^^. Therefore the integral without 
the prefactor gives {A/fi^), and this is arbitrarily large. 
The product of these two factors gives Wtot = 1, exact 
dragging of inertial axes, as demonstrated in Sec. IV Cl 

The de Sitter universe is the limit {p/p) \ (~1) of 
FRW universes with K = 0,±1. These three cases cor- 
respond to three different slicings of one de Sitter space- 
time geometry or part of it. 

How does Mach's principle work arbitrarily close to the 
de Sitter limit, i.e. for [{p/p) -I- 1] <C 1 and p/pcrit finite 
and fixed, when the energy currents within the Hubble 
volume can be made arbitrarily small for any observer? 

For K ~Q we have p? = p/^ — 1&itG{p + p). Hence the 
prefactor in W(x) is of order p^ <^1. On the other hand 
the integral for Wtot would be quadratically divergent 
without the cutoff at Xcutoff = P~^- Hence the integral 
without the prefactor is of order l/p^ ^ 1. The two 
effects cancel, and Wtot = 1 as stated in Eq. pT|) . 

For K — —1 the discussion of the de Sitter limit, (p + 
p) — !■ 0, is exactly the same as the discussion given above 
for the Milne limt. 

For K = ^-1 and {p + p) ^ p, i.e. /2 ^ 1, one obtains 
p ±2i{l - pV8), and W(x) = 8(37r)-i sin"* x, i-e. 
drops out. This gives W(x)dx — 1- 



G. Globally rigid rotations: Zero modes of (A -f 4_R') 

We have seen in Sec. II Bl that for vorticity perturba- 
tions the fixed-time slices St are uniquely determined, 
and that each slice has the unperturbed 3-geometry of 
{E'^, S^, H^}. For the unperturbed 3-geometries we have 
chosen unperturbed 3-coordinates: spherical coordinates 
for {E^,S^,H^}, resp. Cartesian coordinates for E^. 
This does not yet fix the coordinatization (gauge) com- 



18 



pletely, because there still exist residual gauge transfor- 
mations by time- dependent globally rigid rotations of the 
spherical coordinates for 3-space. In this subsection we 
show that these residual gauge transformations are zero 
modes of (A + AK), where — 1. 

Since this paper on Mach's principle focusses on ro- 
tational motion of fiducial axes, not on linear accelera- 
tion of frames, we apply a change of the spatially global 
(rigid) reference frame such that the 3-coordinate veloc- 
ity dx'^/dt of the reference FIDO at the spatial origin, 
X = 0, is left unchanged. The angular velocity of the 
new reference-FIDO axes at the spatial origin relative to 
the old ones is denoted by [—fl*{x = 0)]. 

In the new coordinates the matter 3-velocity field is 
changed by the additional field Wmattorj which is equal to 
the negative of the additional shift vector field /3* and also 
equal to the negative of the additional gravitomagnetic 
vector potential A*, 

Matter = = f^* (0)60. (85) 

The momentum constraint, Eq. ([5]), for the additional 
fields Vmattor = ' with the definition of fl^ for the 
coefficient in the source term, Eq. ([BS]). and with the 
relation {fi^ - M^) = 4if, Eq. dTD]), gives 

(A 4- 4K)A* = 0. (86) 

i.e. the additional field A* arising from a globally rigid 
rotation is a zero mode of the operator (A -I- AK). 

The curl of gives the additional gravitomagnetic 
field B*, 

Bl ^ -2r2*(0)cos6i 

m = +2r!'^(O)sin0 4^ 
" dx 

B; = 0. (87) 

For K = Q (Euclidean 3-space) the additional _Bg -field 
and hence the f2*yj.Q-field are homogeneous fields. But 

for K = ±1 the additional i3g*-field is not a homogeneous 
field, its magnitude varies as a function of x- 

The additional fields Wmattori j B* are equal to the 
fields, which arise in an unperturbed FRW universe, if the 
coordinate system is in rigid rotation relative to FRW 
coordinates, and if the FIDOs and their LONBs are 
adapted to the rotating coordinates. — But note that 
the constraint for a given vanishing source field = 
has the unique solution Ag = 0. 

H. Six fundamental tests 
for the principle formulated by Mach 

The following six tests cover different aspects, inti- 
mately related, of the principle formulated by Mach: 



1. Exact dragging of local inertial axes by a weighted 
average of cosmological energy currents. The 
weight function must be normalized to unity. 

2. "Totally determined" : The time-evolution of local 
inertial axes (relative to any given set of axes) is 
totally determined by cosmological energy currents. 

3. No absolute element as an input. The time- 
evolution of local inertial axes must be an output 
determined by cosmic energy currents, not an in- 
put. 

4. "Solutions Frame-Invariant": The solution for 
rigyro in terms of cosmological energy currents must 
be form-invariant under transformations to glob- 
ally rotating frames, i.e. frame-invariant. 

5. No boundary conditions at spatial infinity in open 
cosmologies are needed to totally determine local 
non-rotating frames. 

6. No reference-frame condition along one world line 
is needed to totally determine local non-rotating 
frames; In open universes such a reference-frame 
condition along a world line at spatial infinity is 
equivalent to boundary conditions. In closed uni- 
verses such reference-frame conditions along one 
world line play the analogous role to boundary con- 
ditions for open universes. 

Tests (l)-(4) represent different aspects of the postulate 
as formulated by Mach. — Test (3) is the postulate of the 
"relativists" (Huyghens, Bishop Berkeley, Leibniz, Mach, 
and others) that only relative motion of physical objects 
is relevant. This was the starting point of Mach. — Test 
(4) is a necessary condition, formulated by Mach in his 
earliest writings [3] quoted in Sec. II Al — Test (1) is an ex- 
plicit formulation of Mach's Principle, given in Sec. II Al 
while Tests (2)- (6) are less explicit about the implemen- 
tation. — Tests (1) and (2) are unambiguous for linear 
perturbation theory. In the nonlinear case, there is now 
agreement that gravitational waves make contributions. 
Hence Einstein's formulation of Mach's principle [s^ . 
"the (metric tensor) field g^i/ is entirely determined ... 
by the energy tensor of matter" is untenable. But the 
definition of "energy currents" ( = "momentum densi- 
ties") and "angular momentum densities" of the gravi- 
tational field has been controversial. Frame-dragging ef- 
fects by "the angular momentum of gravitational waves" 
has been investigated in special, non-cosmological mod- 
els in [H, HI]. — Tests (3) and (4) are meaningful in the 
presence of gravitational waves, even if the concept of en- 
ergy currents of the gravitational field is not available in 
a precise way. — Tests (5) and (6) are intimately related. 

The kinetic energy currents, i.e. the LONB compo- 
nents can be measured without prior knowledge of the 
local non-rotating frames. Hence no absolute element is 
needed in the input for our approach. — In contrast the 
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canonical energy currents, i.e. the coordinate-basis com- 
ponents T\ cannot be measured without the prior knowl- 
edge of the local non-rotating frames all over the universe. 
Hence using coordinate-basis components T*^. means that 
one needs an absolute element as an in put , e.g . in the 
approach emphasized and advocated by [ij, [I^ . 

The six tests listed above are not fulfilled explicitly 
(manifestly), if one uses the momentum constraint for 
the coordinate-basis components , as done and advo- 
cated e.g. in [H, [il, [3. These papers give no proof 
that any one of the six tests is satisfied. 

Einstein's theory of General Relativity applied to the 
solar system fails tests (l)-(5). Hence it is not a "theory 
of relativity" in the sense of the relativists before Ein- 
stein, as has been recognized for a long time. 

In this paper we have shown that the six fundamental 
tests listed above are manifestly (explicitly) fulfilled for 
linear perturbations of FRW universes in Cosmological 
General Relativity. 
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APPENDIX: THE DE RHAM-HODGE 
LAPLACIAN FOR VECTOR FIELDS IN 
RIEMANNIAN 3-SPACES: CONNECTING THE 
CALCULUS OF DIFFERENTIAL FORMS AND 
THE CURL NOTATION 

The material in this appendix (in addition to Sec. II E[) 
is essential for cosmological gravitomagnetism, because: 

(1) For vector perturbations the de Rham-Hodge Lapla- 
cian ("the" Laplacian in the mathematics literature) 
has many crucial advantages (listed below) over V^, the 
"rough Laplacian" , which, unfortunately, has been used 
in the entire literature on cosmology. The calculus of dif- 
ferential forms is needed for the dc Rham-Hodge Lapla- 
cian. In addition, this calculus makes computations sim- 
ple, since no Christoffel symbols are required for compu- 
tations. Unfortunately this calculus has not been used in 
the literature on cosmological vector perturbations. 

(2) On the other hand the elementary notation for vec- 
tor calculus in Riemannian 3-spaces with curl and div 
is indispensable for physical insight in gravitomagnetism, 
i.e. for seeing directly that the mathematical structure of 
the momentum constraint for cosmological gravitomag- 
netism on FRW backgrounds with K = (±1,0) is ex- 
plicitly form-identical (without curvature terms) to the 
familiar structure of Ampere magnetism in Euclidean 3- 
space (except for the Yukawa cutoff at the _ff-dot radius) . 

(3) Unfortunately almost all books on differential ge- 
ometry and on General Relativity, and all research papers 



on cosmological perturbation theory do not treat the top- 
ics of the de Rham-Hodge Laplacian and the Weitzenbock 
formula. The latter gives the difference (A— V^)A, which 
is needed when comparing with the literature on cos- 
mological vector perturbations. In books about differ- 
ential geometry these topics are treated in [2^ — 
But, as far as we know, no book on General Relativ- 
ity or differential geometry and no paper on cosmologi- 
cal perturbations gives the explicit equations which con- 
nect the notation of differential forms with the elemen- 
tary notation for vector calculus in Riemannian 3-spaces: 
(curla)p — (★da)^, Eq. (|A.2[) . and for vorticity fields 
(Aa)^ = — (curl curla)^ — —{*d * da)^, Eqs. (|A.10[ 
and lA.lIj) . Ref. [1^ does not give such equations, and 
Ref. [58] instead gives a "rough, symbolic identifica- 
tion" , which cannot be used literally in the equation 
Ad — —curl curia, see the paragraph following the iden- 
tity (|X2)) . 

Our aim in this Appendix is to give the derivations of 
the tools and results needed in cosmological gravitomag- 
netism. We present this for cosmologists with little or no 
experience in the calculus of differential forms. 

The identities for vector calculus in Riemann 3-spaces 
at the level of first covariant derivatives V are indistin- 
guishable from identities for vector calculus in Euclidean 
3-space, given e.g. in [s^], because curvature effects can- 
not appear at the level of first derivatives. — At the level 
of second covariant derivatives V, the standard identi- 
ties of vector calculus in Euclidean 3-space, e.g. in [STj . 
remain true in Riemannian 3-space, if and only if one 
uses the de Rham-Hodge Laplacian A as opposed to the 
'rough Laplacian' V^. 

The reasons which single out the de Rham-Hodge 
Laplacian in the vector sector on curved 3-space, and 
which give its conceptual motivation in physics: 

1. If for a vector field all types of sources (div and 
curl) are zero, then the de Rham-Hodge Laplacian 
of this vector field is also zero, i.e. the vector field 
is harmonic; Sec. II El 

2. The de Rham-Hodge Laplacian commutes with 
curl, div, grad; Sec. II El 

3. The identities of vector calculus in Euclidean 3- 
space [13] and the important first identity of Green 
generalized to vorticity fields, Eq. l|A.9p . remain 
true in Riemannian 3-spaces, if and only if one uses 
the de Rham-Hodge Laplacian. 

4. For electromagnetism in curved space-time the 
equivalence principle forbids curvature terms, if 
and only if the de Rham-Hodge Laplacian resp 
d'Alembertian is used; Sec. II El 

5. The action principle for Ampere magnetism in Rie- 
mannian 3-space, Eq. (jA.6[) . which is bilinear in 
first derivatives, produces the Ampere equation 
with the de Rham-Hodge Laplacian and without 
curvature terms, Eq. (|A.10|1 . 
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6. The momentum constraint for cosmofogical grav- 
itomagnetism has identical mathematical structure 
to Ampere magnetism (except for the Yukawa cut- 
off), if and only if the de Rham-Hodge Laplacian is 
used; Eq. (O. 

7. The Hodge decomposition for cfosed Riemannian 3- 
spaces, speciahzed to 1-forms, states that any vec- 
tor field can be uniquely decomposed into a sum of 
a gradient field plus a curl field plus a harmonic 
field. This theorem is valid, if and only if the 
de Rham-Hodge Laplacian is used to define har- 
monic fields. 

8. No Christoffel symbols resp connection coeffiecients 
are needed to compute the de Rham-Hodge Lapla- 
cian; Eqs. (lA31lA31[Al4llA:75|l . 

Every one of these properties does not hold for the "rough 
Laplacian" V^. 

1. Differential forms and curl, grad, div 

One can represent a vector field v in Riemannian 3- 
space, a geometric object, by its Local Ortho-Normal 
Basis (LONB) components, denoted by a hat over the 
index, w^, or by its contravariant components in a coor- 
dinate basis w^, or by its associated 1-form v, resp. by 
its 1-form components (covariant components) v\. The 
calculus of differential forms is fundamentally tied to the 
representation in some unspecified coordinate basis and 
with covariant components. 

For the vector product ( a x 6 ) the two notations are 
connected by the equation 

{axb)x= gx^e'""'a^,b^ = {ir[a/\b])x, (A.l) 

where a and b are the 1-forms (covariant vectors) as- 
sociated with a and 6, i.e. (a)^ — (a)^ — a^. A 
tilde designates p-forms (totally antisymmetric covari- 
ant tensors of rank p). — The first step is the wedge 
product A (exterior product, antisymmetric product) of 
two 1-forms, which gives a 2-form c with components 
c^iv = [o,/\b]^i, = a^bv — ai,b^. — The second step is tak- 
ing the Hodge dual (Hodge star operator) of the resulting 
2-form. All our bases will have right-handed orientation. 
The primary definition for the Hodge dual is given in 
any LONB: One contracts with the totally antisymmet- 
ric e-tensor, the Levi-Civita tensor, in a LONB, where 
ej23 = +1- Explicitely: the Hodge dual of of a p-form c 
is given by (*5)j,,.j_^ = (VpO en,..,j„_,j\,..,j, %,..3^, 
where n = 3 is the dimension of our Ricmann space. Ex- 
ample: ft.'-^-' = written in LONB components gives 
/ij = &23- It follows that = 1 for p- forms in a Rieman- 
nian 3-space. — The volume element spanned by a LONB 
with postive orientation is equal to -1-1: ?;(ej, 6*2,63) = 
e(ej, 6*2, 63) = ej23 = +1- The LONB components of the 
e-tensor are invariant under proper rotations (as are the 



LONB components of the metric tensor, g~^j = Sij). — 
But in the calculus of differential forms the Hodge dual 
involves the totally antisymmetric £-tensor in a coordi- 
nate basis, ex^iu = £xp.o V5 and e^^"" = ej^^- (^s)"^- The 
volume of a parallelepiped spanned by a triple of vectors 
is v{X,Y,Z) — e\fi,^X^Y^ Z" , hence the volume 3-form 

V has components wa^i^ = Sx^v — For the Hodge dual 
which takes a p-form into a (3 — p)-form, one can either 
first raise the indices of the original p-form and then act 
with EXfiu, or equivalently one can first act with g-^A"^ a,nd 
then lower the indices to obtain the final (3 — p)-form. 

The operation curl takes a vector field a into another 
vector field of opposite parity, ( curia) — ( V x a ), where 

V is the covariant derivative. The formula analogous to 
Eq. (|A.1|) gives, 

(curla)A = (V X a)x = ffA«e"^"5^az. - i*dd)x. (A.2) 

The last and the next-to-last expressions in Equa- 
tions (jA.2p give the simplest method to compute a curl. 
But the first and second expressions are crucial to make 
the precise connection with the notation universally used 
in electrodynamics via an identity in the form of 
an explicit equation. — The first step on the right-hand 
side of Eq. (|A.2[) . the exterior derivative (antisymmet- 
ric derivative) d of a 1-form d gives the 2-form dd with 
components {dd)^y = df^a^, — d^a^. Because the calcu- 
lus of forms is tied to (unspecified) coordinate bases and 
to totally antisymmetric covariant tensors, the covariant 
derivative V^, after antisymmetrization, can be replaced 
by the ordinary partial derivative 9^. This is the first 
great simplification, and it takes place, because the con- 
nection coefficients in the coordinate basis, the Christof- 
fel symbols F'^^^^, are symmetric in the second and third 
indices. — Compared to other identities given for curl in 
the literature [2J|, the identity (curla)^ = {-kdd)x has the 
important advantage that it directly transcribes the com- 
putational prescription of elementary vector calculus for 
Cartesian coordinates, (curla)^ = SijudjOk, to Rieman- 
nian 3-spaces (or to curvilinear coordinates in Euclidean 
3-space). — Some authors, e.g. [ss], make a 'dictionary' 
with the 'rough, symbolic identification' curlA <^ da^ , 
where the 1-form is the covariant expression for the 
vector A. This identification cannot be used literally in 
the operation (curl curl A), because the Hodge star op- 
erator is missing in front of da^ in their rough, symbolic 
identification. The operation (curl curl) shows that the 
output of the operation curl cannot be the 2-form da^ , 
the output must be a vector (represented e.g. by a 1- 
form) as in our identity (|A.2[) . which directly produces 
the identity (curl curia )a = d* dd)x- 

The operation grad takes a scalar field 4> into \7(j). 

(grad./)) A = iV(b)x = dxq^ = {dcj))x. (A.3) 

The inner product (scalar product) of two vectors can 
be rewritten as an antisymmetric product by first going 
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to the dual of b, 

(A.4) 

where we have used ^Sa/s-yS^''^ = 6^. The wedge product 
of a 1-form d with a 2-form c is given by [d A c]\^i, — 
axCfj^i, + a^Cux + a^cx^ = [c A a\x^v Similarly the inner 
product of two p-forms a, [3 is given by < a, /3 >= *(a A 

The divergence, div a = V • a , can be written as an an- 
tisymmetric derivative in analogy with Eq. (|A.4p . This 
is the second great simplification, because the antisym- 
metrization again eliminates the need for Christoffel sym- 
bols. It provides the easiest way to obtain the explicit 
expression for ( V • a ) at the end of the following sequence 
of equations, 

div d — V ■ d — i< di^ d 

= e'^^-'d^iep^sa') = (1/Vff) (A.5) 

Some well-known textbooks on General Relativity need 
more than a dozen steps to derive the well-known last ex- 
pression in Eq. (jA.sp . because they work with Christoffel 
symbols. — But note that the expression {-kd-kd) is much 
more useful, since it gives e.g. div of Eqs. (|441 [46l) in 
a trivial way. 



2. The Laplacian on vector fields in curved 3-spaces 

The Laplacian A acting on vector fields resp 1-forms 
in curved space is called de Rham Laplacian in js^], and 
[sst writes: "This Laplacian was defined first by Kodaira 
and independently by Bidal and de Rham." The his- 
tory including references is given in ^59J . footnote on 
p. 125. Since other authors use the name Hodge in 
this connection, and in order to be clearly understood, 
we shall call it the de Rham-Hodge Laplacian. The de 
Rham-Hodge Laplacian A must be dis ting uished from 
the "rough Laplacian" = V"Va H^, [IJ. Only when 
applied to a scalar field, one has (A — V^)(/) = 0. — 
For the primary conceptual definition of the Laplacian 
acting on vector fields, we start from an action integral, 
because an action integral only involves first derivatives, 
hence it is independent of curvature effects. For Ampere 
magneto statics in static curved 3-space, the action inte- 
gral (times —1) reduces to the energy functional, 

E[A] = J dv [(curU)^ - An A- J], (A.6) 

where dv = d'^x ^ is the invariant volume element. E\A\ 
is given by the global scalar product of two vector fields, 
(.. , ..) = / dv (.. , ..), where (.. , ..) is the point-wise scalar 
product. 

To carry out the variation of the energy functional 
E\A\ with respect to variations of A, one has to perform 



a partial integration. — To make the surface terms disap- 
pear, we must either have a compact Riemann space, or, 
in a non-compact Riemann space, one of the two forms, 
the variation, must have compact support. — For partial 
integration in curved space one needs the tool of differ- 
ential forms. The adjoint operator to d in the sense of 
the global scalar product for compactly supported forms 
is denoted by d* . Hence {(3P,daP-'^) = {d* (3^ ,aP-^), 
where the superscript p refers to p-forms. Expressing the 
point-wise inner product via a wedge product, Eq. (jA.4|) . 
and using the explicit coordinate-basis expression, partial 
integration gives 

d* aP = {~IY kd* aP. (A.7) 

The superscript * denotes the adjoint, the non- 
superscripted ★ denotes the Hodge dual. — The codijfer- 
ential operator S in this paper is defined with the opposite 
sign from the adjoint, 

S = "d* ^ (5a = diva, (A. 8) 

where a is a 1-form. Some books define d with the op- 
posite sign from our sign, 6 = d* . — The exterior deriva- 
tive d takes a p-form into a, {p + l)-form, the codifferen- 
tial operator S takes a p-form into a (p — l)-form, and 
dd = 0,SS^ 0. 

Green's first identity generalized to vector fields with 
zero divergence on Riemannian 3-spaces (partial integra- 
tion) reads 

J dv (curl a) • (curl (3) ^ ~ J dv a ■ A(3 (A.9) 

in the absence of surface terms. This identity is form- 
identical with the identity in Euclidean 3-spaces, and it 
can be used to give the primary operational definition of 
the Laplace operator (de Rham-Hodge Laplacian) on vec- 
tor fields with zero divergence in Riemannian 3-spaces. — 
Equivalently, in our physics situation the primary con- 
cept of the Laplacian is given by starting from the energy 
functional Eq. (|A.6P : The factor multiplying the varia- 
tion of A in the integrand of Eq. (jA.6P gives Ampere's 
law for A, i.e. AA in the vorticity sector, 

AA = -curl curU = -AnJq (A. 10) 

{AA)^^-{*d*dd)^. (A.U) 

Equation (jA.lip with the operator {-kd-kd) gives the 
simplest computational method in Riemannian 3-spaces 
(and for curvilinear coordinates in Euclidean 3-space). 
On the other hand, Eq. (|A.10p with the operator (curl 
curl) in Riemannian 3-spaces is needed to make the pre- 
cise connection using the notation universally adopted in 
electrodynamics in Minkowski space in a (3-t-l)-split [43 |. 

In the general case of vector fields with div a ^ and 
curl a ^ the requirements for the Laplace operator A 
of de Rham and Hodge acting on vector fields are: 

(1) 2nd order elliptic differential operator, 
self-adjoint, (A.12) 

(2) curl a = anrf div a = => Aa = 0. (A.13) 
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Reversing the arrow in the statement (jA.lSp would be in- 
correct for open Riemannian 3-spaces (open universes) in 
contrast to statements in some textbooks. This is evident 
from simple counter-examples: For a FRW background 
with K — a, homogeneous i?-field gives AA — 0, but 
curl A ^ 0. For a FRW background with K — —1 the 
analogous field is a poloidal vorticity field B with £ = 1 
and q = 0, for which again A A = 0, but curl A ^ 0. — 
From Eqs. (|A.12| IA.13"1) it follows that 

Aa = —curl curl a + grad diva. (A. 14) 
Aa — d-k d + d-k d-k) d 

= {Sd + dS)d 

= -{d + d*fh. (A.15) 

Note that {dd*) is self-adjoint, and [d+d*f = {dd*+d*d) 
is self-adjoint and positive semi-definite. 

The "rough Laplacian" = g^^V i, does not satisfy 
the requirement {curl a = and div a = 0} Aa = in 
curved space, and it does not satisfy Green's first identity 
Eq. ((X9ll . 

3. The Weitzenbock Formula 

The Weitzenbock formula [2^, [g^] for vector fields 
[2^ 1 and in 3-dimcnsional Riemannian space is de- 



rived by first using (curl curl A)a = — (V^A)q -t- 
(grad d\vA)a + [(V/jVa - VaV/3)^]''. For the last term 
one uses the Ricci identity, 

[(V^V. - V.V^)^]" = R%^A^ (A.16) 

and the contracted Ricci identity, 

[(V^V, - V,V^)v4]^ = R.pA^, (A.17) 

where R^p^^, is the Riemann tensor and Ri,f} the Ricci 
tensor. The sign conventions for the Riemann and Ricci 
tensors are those of Misner, Thorne, and Wheeler [20j . 
This gives the difference between the true Laplacian A 
(de Rham-Hodge Laplacian) and the "rough Laplacian" 
V"Vq,, i.e. the Weitzenbock formula, 

{AA - V2i)„ = RjAf,. (A.18) 

For a FRW universe with K ^ ±1 and a curvature scale 
Qc, the Ricci tensor of 3-space is R£ — {2K/a1) 6£ , and 
the Weitzenbock formula gives 

{A-V^)A^-{2K/al)A. (A.19) 
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